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GENERATION Or KONG DISPERSE DROPLETS BY SPONTANEOUS 
CONDENSATION OF FLOW IN NOZZLES 


Abstract 


Submicrcn size mcnodisperse particles are of interest in many 
industrial and scientific applications. These include the 
manufacture cf ceramic parts using fine ceramic particles, the 
production cf thin films by deposition of icnized c lus t- cy r . , 
mcnodisperse seed particles for laser anemone t ry , and the study 
cf size dependence cf cluster chemical and physical properties. 
An inexpensive and relatively easy way to generate such particles 
is by utilizing the phenomenon of spontaneous condensation. The 
phenomenon occurs when the vapor or a mixture of a vapor and a 
noncondensing gas is expanded at a hich expansion rate. The 
saturation line is crossed with the supercooled vapor behaving 
like a gas, until all of a sudden at so called the Wilson point, 
condensation occurs, resulting in a large number of relatively 
mcnodisperse droplets. The droplet size is a function of the 
expansion rate, inlet conditions, mass fraction of vapor, gas 
properties, etc. 
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PRECEDING BLANK NOT FILMED 



In this dissertation, spontaneous condensation of (a) steam and 
(b) water vapor and air mixture in one dimensional nozzle has 
been modeled and the resulting equations solved numerically. The 
classical (Frenkel and Volmer) and Deich's nucleaticn theories 
were used. Unlike classical theory, Deich's theory also accounts 
for the critical droplet size dependence upon the expansion rate 
in the nozzle. The droplet size distribution at the exit of 
various one dimensional nozzles and the flow characteristics such 
as pressure ratio, mean droplet radius, vapor and droplet 
temperatures , r.ucleation flux, supercooling, wetness etc., along 
the axial distance were obtained. The numerical results compared 
very well with the available experimental data. Beth nucleaticn 
theories predicted the pressure ratio properly but the classical 
theory tended to overestimate the droplet sizes. Additionally, 
since the classical theory cannot be used for obtaining a 
correlation for surface tension ratio, Deich's theory is 
preferred for analyzing and designing new nozzles. 

The effect of inlet conditions, nozzle expansion rates and vapor 
mass fractions cn droplet mean radius, droplet size distribution 
and pressure ratio were examined. By manipulating the various 
parameters, different mean radius and pressure ratio can be 
obtained. However, this only has minor impact on the relatively 
monodisperse droplet size distribution. 
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CHAPTER 1 


INTRODUCTION 


Submicron size monodisperse particles are of interest in many 
industrial, as well as scientific applications. These include the 
manufacture of ceramic parts using ceramic particles, the 
production of thin metal, semiconductor and organic films by 
deposition of ionized clusters, and the study of the size of 
dependence of cluster chemical and physical properties [1.01]. 
Monodisperse seed particles are also of great importance in laser 
anemometry. The ultimate accuracy of particle scattering based 
laser velocimetry measurement depends upon the ability of the 
light scattering seed particles to faithfully track the local 
flow field. For achieving maximum measurement efficiency, the 
control of seed particle size, size distribution, shape and a 
specific material density is desired [1.02]. Edwards [1.03] 
suggests that all possible bias problems in laser anemometry 
would be easier to be dealt with if there existed an inexpensive 
seeder that was capable of generating monodisperse seed particles 
fast enough to provide acceptable levels of data density. The 
size ranges of 0.1 to 1 /am are of great interest in seeding high 
speed flows such as high subsonic, transonic, supersonic and 
hypersonic flows. 
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The phenomenon of spontaneous condensation, also called 
homogeneous condensation, can be utilized for generating 
monodisperse particles. The phenomenon occurs when the vapor or 
the mixture of vapor and gas is expanded at a high expansion rate 
through a nozzle, such that the saturation line on the Mollier 
chart is crossed with the supercooled vapor behaving like a gas 
until all of a sudden at a certain point called the Wilson point, 
condensation occurs, and the condensate separates out, resulting 
in a large number of nearly monodisperse droplets. The droplet 
size is a function of the expansion rate, nozzle geometry, inlet 
total conditions, mass fraction of vapor in case of vapor gas 
mixture, gas properties, etc. 

The phenomenon of spontaneous condensation is also of much 
interest in the problem of wetness loss in the low pressure as 
well as high pressure stages [1.04], [1.05], [1.05] of steam 
turbines. The latter case was stimulated by the advent of the 
light-water-cooled reactors which supplies dry saturated steam at 
inlet to the high pressure cylinder of the turbine. The same 
phenomenon was found important in the supersonic and hypersonic 
wind tunnels [1.07], [1.03], [1.09] and free-jet expansions 
[1.10]. Condensation is usually undesirable in these 
applications, and therefore a bezter understanding and prediction 
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of the phenomenon are crucial to the control or even elimination 
of the condensation. 

Experimental and theoretical studies of expansion of steam 
through supersonic nozzles has been conducted by various 
investigators [1.11], [1.12], [1.13], [1.14], [1.15], [1.16]. 
Experiments conducted by Gyarmathy and Meyer [1.13] show a strong 
relationship between the rate of expansion in the nozzle and the 
spontaneous condensation process. The theoretical investigation 
of steam nozzle flows [1.13], [1.14], [1.16], [1.17] have mostly 
been conducted from the point of view of obtaining the pressure 
rise, moisture levels and the decrease in the nozzle efficiency 
associated with the spontaneous condensation process. More 
recently, Turner et al. [1.01] have investigated the use of 
subsonic nozzles for the formation of monodisperse particles in a 
system of interest for industrial powder production. A mixture of 
nitrogen and aluminum sec-butoxide vapor was considered. 

While several theories have been proposed for obtaining the 
nucleation rate during spontaneous condensation, the two main 
theories of interest are the classical theory (known as the 
Frenkel and Volmer theory) and the theory developed by Deich et 
al . The most widely adopted classical nucleation theory is 
discussed by many authors [1.18], [1.06], [1.19], [1.20]. Deich's 
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nucleation theory is reviewed by Moore [1.16] and Kadambi and 
Whir low [1.14]. The derivation of the classical theory is based 
upon the kinetic theory [1.18], The nucleation rate can be 


expressed as a 

function of 

thermodynamic properties 

and 

the 

surface tension 

coef f icient . 

Deich's theory accounts 

for 

two 


components of the formation of nuclei. The first one arises from 
the so-called " heterophase fluctuations", which means that 
statistically a number of molecules will clump together and form 
a droplet size equal to or greater than the critical size 
necessary for survival and subsequent growth. This component 
corresponds to the classical theory but considers more realistic 
conditions about the steadiness of the droplet generation and is 
considerably smaller in magnitude than the classical theory. The 
second component arises from the fact that the critical size 
depends upon the vapor condition in a fashion such that during 
expansion in the nozzle, the critical droplet size decreases. As 
a result, some droplets formed at subcritical sizes are 
transported, before they can evaporate, to a downstream region 
where their sizes become super critical. The final nucleation 
rate is the addition of the two components. It can be shown that 
the second component is usually of order of magnitude much higher 
than the first one and is directly proportioned to the expansion 


rate in the nozzle. 
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One of the major drawback of the classical theory is that it does 
not directly account for the effect of the rate of expansion on 
the nucleation rate, while Deich's theory does. According to 
Gyarmathy [1.13], the experimental observation shows that 

expansion rate has strong effects on the nucleation rate. It is 
also noted that all nucleation theories cannot achieve 

satisfactory agreement with experimental data without using any 
correction factors. One of the unknown parameters in the 
nucleation theory equations, which also has substantial effects 
on the nucleation rate, is the surface tension for small 
droplets. By using surface tension ratio as correcting factor, 
Moore [1.16] showed that the classical theory cannot obtain 
agreement with the experimental data of pressure ratio and mean 
droplet radius simultaneously. The other drawback of the 
classical theory is that it can not establish usable correlations 
between the flow conditions and the correcting factor. Such 
correlation is especially desirable if the method has to be 
utilized for predicting various new flow conditions. Salcanov 
[1.17] has used surface tension for flat surface in his study and 
uses a correction coefficient in the exponential term for the 
nucleation rate in the classical theory equations. He varied the 
value of the coefficient to obtain agreement v/ith experimental 
pressure ratio and obtained a correlation between the correction 
coefficient and the inlet total pressure. However, only three 
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data points were considered in a pressure range < 40 psia and 
there was considerable scatter. 


1 . 1 Obj ect ives 


The objective of this work was to develop a theoret ical /numer ical 


model of 

the flow 

of vapor and vapor/ inert 

gas 

mixture 

in 

a 

supersonic 

nozzle 

with 

spontaneous 

condensation, 

capabl 

0 

of 

predict ing 

droplet 

size 

distribution , 

The 

mode 1 

wou Id 

be 

a 


useful tool to design nozzles to obtain desired micron sized 

droplets. Steam was selected as the vapor because of 

availability of experimental data to verify the model. The 

classical and the Deich's nucleation theories were used in this 


study to provide a comparison between the two. 



CHAPTER 2 DESCRIPTION OF SPONTANEOUS CONDENSATION PHENOMENA 

2.1 Introduction 

To generate relatively monodisperse particles, a spontaneous 
condensation based monodisperse particle generator is proposed. 
The analysis will consist of the expansion process of both the 
supercooled vapor phase and condensed liquid phase. The variables 
which define the character ist ics of the two phases are related by 
the conservation laws of mass, momentum and energy of the whole 
system as well as the “life history 1 * of the condensate: when are 
the droplets born? what are their sizes and number when first 
born? how will they grow? Before we present all the governing 
equations and calculation procedures, it is beneficial to define 
some of the fundamental concepts and outline the specific problem 
we will be dealing with in the later analysis. 


2.2 Homogeneous Condensation In A Supersonic Nozzle Flow 

Let us consider a vapor expanded from a superheated state such 
that it crosses the saturation line. If the expansion is 
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relatively fast, the vapor will not have enough time to adjust 
itself immediately for the condensation. Rather, it will remain 
as vapor even after crossing the saturation line. This vapor 
phase which is in a meta-stable equilibrium condition and behaves 
like superheated vapor is referred to as “supercooled vapor*. 
Very often, the name "supersaturated” or "subcooled " is also 
used . 

When a vapor or a mixture of vapor and gas is expanded in a 
nozzle in a relatively high expansion rate, condensation may be 
delayed with respect to the equilibrium states (as described in 
the previous paragraph, no condensation occurs even after 
crossing the saturation line) . The supercooled vapor continues to 
expand just like a superheated vapor. At certain point called 
Wilson point where the supercooled vapor can no longer hold its 
meta stable equilibrium condition, further expansion causes a 
sudden collapse of the supersaturation with the formation of 
nearly monodisperse droplets appearing in the form of a fog. This 
process is called spontaneous condensation or homogeneous 
condensation which takes place in absence of any foreign nuclei. 
This phenomenon is also observed in supersonic wind tunnels (when 
humidity is high and the moisture separates out all of a sudden 
in the form of a fog) , wet steam nozzles and steam turbine blade 
passages. The release of latent heat during condensation results 
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in a diabatic flow with an increase in pressure. 


The phenomenon for steam or any other vapor flow in a condensing 

Laval nozzle is shown in Figure 2.1. The enthalpy-entropy diagram 

shows the isentrcpic expansion vapor from initial stagnation 

conditions to crossing the saturation line and staying in a 

supersaturated state until at pressure p when suddenly 

w 

spontaneous condensation occurs. The release of latent heat 
causing an increase in temperature and pressure as well as 
entropy. This pressure rise is often misinterpreted as the 
condensation shock. As will be shown later, the steepness of the 
pressure rise depends mainly on the expansion rate as well as the 


staanat ion 

conditions of the 

vapor 

prior to 

expansion. At 

low 

expansion 

rate the pressure 

rise 

is not 

appreciable . 

The 

supercooling decreases as the 

vapor 

temperature, pressure 

and 


entropy increase. Consequently, the supercooled vapor becomes 
closer and closer to its equilibrium state. The spontaneous 
condensation process takes place in a very short time. 
Eventually, at the end of spontaneous condensation shown by 
pressure in Figure 2.1, the supercooled vapor will almost 

achieve its equilibrium state and further expansion would take 
place almost like an equilibrium isentrcpic expansion again. 
There will still be a very small amount of supercooling at the 


end of spontaneous condensation, and the supercooled vapor will 
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gradually recover back to the equilibrium state as the expansion 
continues or the adjustment would occur in the form of either 
weak oblique shock cr expansion wave at the exit of the nozzle 
depending on whether the atmosphere pressure is smaller or larger 
than the exit pressure. Since the degree of supercooling at this 
stage is so small (usually only a few degree F) , the deviation 
from the equilibrium case is usually negligible. Besides, the 
complex non-equilibrium phenomena occurs at the exit of the 
nozzle is not our iriain interests thus will not be oiscussed in 
this work. 

A lot of information which would serve for different applications 
can be obtained from the simulation of the process described 
above. For example, the information about velocity and Mach 
number variation, as well as the Wilson point and the pressure 
distribution would be useful in wind tunnel application; the 
wetness level would be useful in examining the erosion on turbine 
blades; the entropy increase would be valuable in calculating the 
loss in total pressure or the efficiency on turbine application; 
etc. Our main interest is to obtain the mean droplet size and 
size distribution and proper correlation between the crucial 
parameter and the inlet condition so that one can design the 
nozzle to give us the desired size and distribution. 
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2.3 Fundamental Definitions 

We have illustrated the physical phencmincn and the model we will 
be analyzing in the last setion. Before we can go on to the 
mathematical simulation equations, we need to quantitatively 
define the degree of supersaturat ion and some of the important 
physical qantities. 


2.3.1 supercooling 

3y assuming isentropic expansion and treating the vapor phase as 
an ideal gas with constant specific heat at constant pressure c 

pg 

and constant specific heat ratio y, we can greatly simplify our 
analysis without losing the essence of the phenomena. As shown in 
Figure 2.2, suppose the vapor is expanded isentropical ly from the 
superheated state noted as point T, across the saturation line at 
point A to a certain supercooled state 3 on the Mollier chart 
(enthalpy-entropy chart) . The location of point 3 on the Mollier 
chart can be determined by the enthalpy drop Ah between the two 


states A and B 
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By shifting the expansion line to the right and left, we can 

eventually construct the isobar passing point B which intersects 

the saturation line at point C as shown in Figure 2.2. For ideal 

gas, the isotherms coincide with the constant enthalpy line. 

Since the supercooled vapor is considered to behave just like the 

superheated vapor, the isotherm crossing point B will be the 

extension of the isotherm in the superheated region. This 

isotherm shown in Figure 2.2 has the intersection with saturation 

line at point D. Therefore, the saturation temperature 

corresponding to the supercooled vapor pressure p will be T (p ) 

B S 3 

= T^ . From Figure 2.2, we can see that 

T < T = T (p ) (2.2) 

B C S 3 

The temperature of the supercooled vapor is always lower than the 
saturation temperature corresponding to its own pressure. The 
supercooling AT is defined to be the temperature difference 
between the two: 

AT = T (p ) - T 
s 3 9 


(2.3) 
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Supercooling represents the deviation of the state of the 
supercooled vapor from the equilibrium condition. The higher the 
supercooling, the larger the extent of the non-equilibrium. 


2.3.2 supersaturation 

From Figure 2.2, we can also see that the saturation pressure 

corresponding to temperature is the same as the pressure at 

point D, p (T ) = p . The isobar passing point D lies below the 

s B D 

isobar passing point C which is corresponding to p_ , therefore, 

B 

we have 

p>p=p(T) (2.4) 

B D s 3 


The pressure of the supercooled vapor is always higher than the 
saturation pressure corresponding to its own temperature. The 
supersaturation s is defined as the ratio: 


s 


P (T ) 
s B 


(2.5) 


The natural logarithm of the supersaturation A is usually more 


useful 
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A 


ln( 


p (T ) ) 
s 3 


2.3.3 critical radius 


The stability of a spherical droplet surrounded by the vapor 
phase was studied by Thompson (1870), Helmholtz (1335) and Gibbs 
(1373) [1.13]. Though their work proceeded largely independently, 
by considering the chemical potential and mechanical equilibrium, 
they found the radius of a droplet in equilibrium with its vapor 
can be expressed as 


2 cr 

cr p RT A 
t g 


(2.7) 


An alternative expression shown by Kirillov et al. [1.06] is 


r 

cr 



(2.3) 


The detail derivation can be found from Kirillov et al. [1.06], 


and will not be presented here. 
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The significance of the critical radius is that it defines the 
minimum size of the droplet able to survive under the given vapor 
conditions. If the droplet radius is below the critical radius, 
the chemical potential of the liquid phase increases, and the 
liquid phase tends to convert to vapor, i.e. the droplet will 
vaporize. On the other hand, if the droplet radius is beyond 
critical, the chemical potential of the vapor phase increases, 
and the vapor tends to convert to the liquid phase, i.e. the 
droplet will grow. If the droplet radius is equal to the critical 
radius, its size will remain unchanged provided the pressure and 
temperature of the vapor phase remain the same. Thus, only those 
droplets which have reached the critical size can serve as the 
condensation nuclei. As shown by Kirillov et al . [1.06], the 
critical radius is primarily affected by the supercooling AT (or 
A which is closely related to AT) and is insensitive to the 


change of pressure. 



CHAPTER 3 


NUCLEATION THEORIES 


Though it is net the purpose of this work to review the 
development and status of the spontaneous nucleation theory in 
detail, it is important to consider the basic concepts and do a 
comparison of the two theories used in this work. The detailed 
historical evolution and discussion can be found in various books 
and articles [3.01], [1.18], [1.20], and it is sufficient to say 

that the the theoretical prediction cf the nucleation rate is 
still in doubt by several orders of magnitude. 

In the absence of foreign nuclei or during the rapid expansion 

where the condensate accumulation on the foreign nuclei is small 

enough to be negligible, the condensation will be delayed with 

respect to equilibrium state. As will be discussed later, the 

spontaneous condensation process can be viewed to start from the 

formation of very small droplets of critical size whose typical 

-9 

order of magnitude is around 10 m. These droplets serve as the 
"condensation center", on which the supercooled steam condenses 
thus causing the droplets to grow. The nucleation theories 
predict the rate of formation of such nuclei providing the number 
cf nuclei per unit time per unit volume or mass. 


16 



17 


3.1 Classical Theory 


The most commonly used classical nucleation theory was derived by 
Volmer [3.02], Frenkel [3.03] using a different approach also 
derived the identical equation. The basic idea is to postulate 
[3.04] that the probable number of molecular clusters formed by 
random collisions is proportional to Boltzmann factor, 
exp(-AS/k), where AS is the entropy decrease associated with the 
formation of such cluster. AS can be related to the work 
associated with the formation of the cluster W as AS = W/T. W can 
be further calculated from the sum of different contributions of 
the change of Gibbs free enthalpy AG associated with the 
formation of such cluster. Furthermore, Farkas [3.05] introduced 
a hypothetical quasi-steady process model, in which clusters are 
thought to be continually removed and replaced by single 
molecules, and there exists a net flux of clusters growing into 
critical size. Detailed derivation is described by Wegener 
[1.13], or in the original work by Volmer [3.02] and Frenkel 
[3.03]. The final expression for the nucleation rate is* 


J = Z exp 

o 



1 


(3.1) 


V 


A' 
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where 


Z 

o 


p 


RT p 
g f 



*) 

^2 _ 167lN(7 V 

1 “ 2 3 

3p (RT ) 
f 9 


A 


In 


Ul_1 

[p (T ) J 

v s g y 


With the same basic concepts mentioned above, different authors 
[1.18], [1.04] derived various revised versions of the nucleation 
theories by modifying some of the physical assumptions and 
mathematical simplifications. Nevertheless, except Lothe and 
Pound [3.06] whose theory predicted 10 1? times higher nucleation 
rate, all the other theories result in nucleation rate of similar 
order of magnitude as Volmer and Frenkel's theory. None of these 
theories, including the theory derived by Lothe and Pound, can 
explain the discrepancy between the experimental data and the 
theoretical prediction. 


3.2 Deich's Nucleation Theory 


Deich et al. [3.07] concluded that the classical theory and 


Lothe- Pound theory are both unsound. They proposed a nucleation 
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theory which consists of two components. The first component is 
due to the so called " het erophase fluctuations". Due to the 
fluctuation of the density in the vapor, two or more molecules 
may collide and attach together for a short period. Some of them 
happen to be droplets of critical size or larger. Two assumptions 
were made during the derivation: 

1. The vapor has an equilibrium nuclei number distribution 
before the expansion. 

2. The expansion is so rapid such that the number 
distribution of the droplets does not have enough time to 
adjust, and thus remains at the initial equilibrium 
distribut ion . 

This component is similar to the classical theory yet considers 
more realistic conditions: a non -equilibrium distribution was 

used instead of the equilibrium one. The second component is 
directly related to the expansion rate. As explained in chapter 
2, the critical cluster size is a function of the state of the 
vapor. Suppose some droplets of subcritical size are formed at 
certain point during the expansion, if the condition of the vapor 
remains the same, these droplets would have been evaporated. 
However, since the supersaturation increases and the pressure 
decreases during the expansion, the local critical cluster size 
may be smaller than the radii of the subcritical droplets formed 


previously. As a result, 


these subcritical droplets are 
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transported, before they can completely evaporate, to a 
downstream region where their size is considered to 
supercritical. These droplets then also serve as the condense! , u 
nuclei and constitute the second component of Deich' s nucleation 
theory. Since there are some mistakes in the original paper by 
Deich [3.07] and the equation presented by Moore [1.16], a 
detailed derivation is described in Appendix A. The final 
equations may be expressed as below: 

Deich ' s_theory_( f irst_compor:ent ) 


J 

l 


3ANpv^RT 
g 

/8rr c 


( T Act /<t - T A 1 


Z l 3T A«r l<r - 2T A 

iQ 

O 

O 

o 

exp 

1 g 0 0 0 

T 

2 T A 

l 0 J 


A o 


(3.2) 


where subscript 0 stands for the condition before expansion. 


Deich ' s_theory_(second_component ) 


J 
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N 

dAT 


r T 1 
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5 


g 

T~ 

9 J 



exp 
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(3.3) 


The total nucleation rate is: 
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J = vj + J 
1 2 


( 3 . * ) 


When expanding from superheated vapor, = 0 , thus the two 

components can be further simplified respectively as: 
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(3.5) 


(3.6) 


In most instances, except in the case of very low expansion rate, 
the magnitude of is several orders greater than J^, [1.14] 
thus the total nucleation rate can be approximated by: 


J = J 


2 


(3.7) 


3.2.1 Relation between the nucleation rate and the 
expansion rate 

One of the important term in the dominant component in Deich's 
nucleation theory is . As we mentioned before, the 
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condensation process is strongly affected by the expansion rate 
p. Therefore, it is interesting to see how these two term could 
be related. For isentropic expansion, the two terms can be 
related [3.03] as: 


dAT 

dt 


T (r - i) 
g r 


pv T 
sg s 



(3.3) 


Detailed derivation is shown in Appendix B. Equation 3.3 shows 
dAT 

chat and therefore nucleation flux J is directly proportional 

dt 

to the expansion rate p. 


3.3 Surface Tension Of The Droplet With Small Radius 

The accuracy which can be achieved in calculating non-equilibrium 
condensing nozzle flow depends ultimately on the accuracy of the 
nucleation theory and the droplet growth. The accuracy of the 
practical computation of the nucleation rate, however, depends on 
the accuracy of the properties of the droplets with small radius. 
The most important property which has very strong effect on the 
nucleation rate is the surface tension or. According to many 
different investigators [3.09], [3.10], [3.11], [3.12], [3.13], 
the value of surface tension for a highly curved surface, cr^ may 
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be considerably different from cr , surface tension for a flat 

00 

surface. Some researchers suggest that <r > <r [3.09], [3.10], 

r oo 

while ethers claim <x < cr [3.11], [3.12], [3.13]. Since no 

r co 

definite conclusion can be made yet, so we adopted a correction 

factor, surface tension ratio tr / a , to obtain agreement 

r oo 

between the calculations for the pressure ratio with the 
available experiment data. 

Although the two nucleation theories are substantially different, 
both need the value of the surface tension for droplets with 
small radius. As will be shown in chapter 6, for both nucleation 
theories, only a small change in the surface tension ratio 
results in substantial change in the pressure ratio, mean droplet 
radius, and droplet size distribution. A higher surface tension 
ratio tends to delay the onset of spontaneous condensation 
(indicated by pressure rise) thus decreasing the Wilson pressure, 
increasing the mean droplet radius and causing a slightly larger 
variation in droplet radius. The results are relatively very 
sensitive to the change in surface tension ratio. It will be 
shown later, although the correlations for surface tension ratio 
of the two nucleation theories have similar trend, the one 
obtained using classical nucleation theory is not an acceptable 
correlation due to relatively large scatter in the correlation 


data . 
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It has to be noted that u is a strong function of temperature 
( cr ( T ) ) . To interpret it properly, the vapor temperature j.‘ , 

*5 

instead of the saturation temperature T (p) , should be used, 

s 

because the droplet has the same temperature as the vapor when it 

is freshly formed. Hence, the surface tension should be cr = 

<r(T ). The equation listed bv Leinhard [3.14] with slight 
g 

modification in the coefficient was- used to account for the 
temperature dependence of cr of water. It may be expressed as: 


o' = 235.3 (1 


1.256 

T ) [1 - 0.625(1 - T }] x 10 

R R 


Nm 


where reduced temperature is defined as 


(3.9) 


T 

T = 2 

r 647.286 


(3.10) 


and T is in K. 
g 



CHAPTER 4 . 


DROPLET GROWTH PHENOMENA 


4.1 Droplet Growth Rate (Heat Transfer Coefficient) 

The preceding chapters treated the birth of the droplets, their 
size and number. The next step of the analysis would be how will 
they grow? We have pointed out in the previous chapter that the 
temperature of the freshly formed droplet is equal to the vapor 
temperature. However, as expansion continues, the vapor 
temperature decreases faster than the temperature of the 
droplets. Thus, a temperature gradient is established between the 
droplet and the surrounding supercooled vapor. Consequently, heat 
transfer occurs between the two phases, which causes the droplet 
to grow. The detailed microscopic mechanism of the condensation 
process of a supercooled vapor is quite involved and beyond the 
scope of this work. However, the process can be simply 
illustrated by the following. 

Let us consider a control volume of the vapor immediately 
adjacent to the droplet in Figure 4.1 with temperature 
distribution T(r) . If the droplet temperature is higher than the 
surrounding supercooled vapor temperature T >T , without Juhav.oh; 
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the detail distribution of T(r), it is reasonable to say that 

T >T(r. ) and T(r ) >T . Hence, heat is transfered from the droplet 
f 1 2 g 

to the adjacent vapor Q { r ^ ) as well as from the adjacent vapor to 
the surrounding supercooled vapor Qtr^). When the adjacent vapor 
rejects more heat than it receives Q(r^)>Q(r ), phase change 
cccurrs and the adjacent vapor condenses. Accordingly, the rate 
of condensation depends on the rate at which the latent heat can 
be carried away from the vapor adjacent to the droplet surface 
into the cooler vapor. If Q(r^)<Q{r^), then the adjacent vapor 
will be heated up and there will have no condensation. 

Similarly, in the case of evaporation, the vapor temperature is 
higher than the droplet surface temperature which is higher than 
the droplet bulk temperature. When the droplet surface receives 
more heat than it rejects, phase change occurrs and the droplet 
surface evaporates. 


The growth rate (positive for condensation, negative for 
evaporation) can be calculated from the energy balance of the 
droplet 


4 3 

pc - Trr 

f 3 


dT 

* 

dtT 


. 2 dr 

* 4 " r p f 3F 


h (T ) - h (T ) 
£ t g g 
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4nr a(T - T ) 

g f 


( 4 . 1 ) 
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The heat capacity term in Equation (4.1) is negligible compared 
with the latent heat released as shown by Gyarmathy [1.04], thus 
the droplet growth equation can be expressed as 

4irr 2 p ^ h (T ) - h (T ) 1 = 4nr 2 a(T - T ) (4.2) 

fdt[gg f f J f 3 

To correctly compute the droplet growth rate, a proper expression 

of the heat transfer coefficient a is important. Since the order 

cf magnitude of the droplet diameter v;hen it is first born 

_ c 

(typically, critical radius of steam r s 10 " m) is smaller than 

c r 

the typical value of mean free path (for steam, 1 s 10 at 
atmospheric pressure), the effect of molecular interaction has to 
be taken into account when evaluating a. As droplet grows, the 
droplet diameter may eventually become greater than the mean free 
path. Hence, the molecular interaction is less important and the 
continuum flow prevails. While the transition from rarefied flow 
to continuum flow is really gradual rather than sharply defined, 
the different flow regimes may be represented by different range 
of Knudsen Number 


Kn 


1 

2r 



(4.3) 


where mean free path 1 can be found from kinetics theory [4.01] 
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1 



(■ 1 . 4 ) 


4.1.1 continuous flow regime 


For Kn < 0.01, the mean free path’ is much smaller than the 
droplet radius, therefore, the equations of continuum mechanics 
apply [4.02]. Since the droplet travels in the same velocity as 
the surrounding vapor, the heat transfer process between the 
droplet and the surrounding vapor can be approximated by free 
convection of a sphere at uniform temperature. The Nusselt number 
can be calculated as [4.03] 


where 
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For most vapors, Pr is of order 1, R 


10 “ 


10 4 J/ (Kg K) , and r 
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- 6 

is no larger than 10 m m most spontaneous condensation 

processes. Suppose the temperature difference and the vapor 

2 . - 2 

temperature are bcth of order 10** K # and Kn is at the limit 10 , 

the order of magnitude of the last term in Equation (4.5) for the 

-4 

most conservative case would be 10 which can be safely 
neglected compared with 2. Thus the Nusselt number will be 


Nu 

m 


a(2r) 

k 
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(4.6) 
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(4.7) 


Equation (4.2) can now be expressed as 


sfr 


p — h (T ) - h (T ) 
f dt g g f f 


1 ^ 


(T - T ) 
f g 


(4.8) 


4.1.2 free molecular flow regime 

Opposite to continuum flow, the other extreme is the free 
molecular flow regime where Kn > 3 . Hill [1.19] has derived the 
energy transfer equation under high Kn condition, Young [1.05] 


rewrites it with an alternative expression as 
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dr 
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f dt 
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g g it 


p(jr + l)c 

— - ■ - - (T - T ) 

/ 2nRT 2 7 

V g 


( A 


From Equation (4.9), the heat transfer coefficient for free 

molecular flow a is 
f m 


p(? + Dc 
P 

/ 27lRT g 2 y. 


(4.10) 


which is practically the same expression presented by Gyarmathy 
in reference [1.09] utilizing the derivation by Gyarmathy [4.04] 
and Kang [4.05] 


a 


f n 



y + 1 3 Pr 
7 4Kn ' 


(4.11) 


However, the expression of the Knudsen number Kn ' in Equation 
(4.11) is misrepresented as 


Kn ' 



4rp 


(4.12) 


in Gyarmathy 7 s derivation. Therefore, symbol Kn' instead of Kn is 


used in Equation (4.11). 
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4.1.3 slip and transition flow regime 


There is no simple analysis which can describe the complex 
molecular interactions occurred in the flow of intermediate 
Knudsen number 0.01 < Kn < 0.1 {slip flow) and 0.1 < Kn < 3 

(transition flow) . A simple analytical model suggested by Young 
[1.05] postulates that free molecular flow is valid within a few 
mean free path of the droplet surface, while continuum theory 
prevails outside this region. The model is illustrated in Figure 
4.2. The interface separating the two flow regimes is at radius 
(r + al) where the suitable value of parameter a is suggested by 
Fuchs [4.06] to be 2. According to this model, Equation (4.8) is 
valid outside the interface, and Equation (4.9) is valid within 
the interface. By replacing T with T^ at the right hand side of 
Equation (4.8), we obtain Equation (4.13). By replacing T^ with 

T at the right hand side of Equation (4.9), we obtain Equation 

i 

(4.14). The two equations are shown as follows: 
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(4.13) 


(4.14) 


Note that the latent heat part remains the same as Equation (4.8) 
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and (4.9), for the transition region in the vicinity of the 
droplet is considered to retain the same condition and propagate 
as the droplet grows. As a result, the net energy change of phase 
transformation is the enthalpy change between the droplet and the 
surrounding supercooled vapor at infinity. Combining Equation. 
(4.13) and (4.14) by eliminating T., we get 

P f ~ |h (T ) - h (T ) 
f at I g g f f 

This equation is identical to the one derived by Young [1.05] 

except the more realistic enthalpy change h (T ) - h (T ) is used 

g g f f 

rather than latent heat, and the expression of Kn is also 
misrepresented as Equation (4.12) in Young's equation. With a = 
0, the equation is identical to that derived by Gyarmathy [1.09] 
except Kn should be replaced by Kn ' . Equation (4.15) net only 
correctly approaches the limit of the two extreme flow regimes, 
but also provides a smooth transition between the two as well. 
Therefore, it is thought that Equation (4.15) can be used for all 
values of Kn. It is also noted by Young [1.05] that Gyarmathy ' s 
equat ion 

a 

c 

“ = 1 + 3 . 18Kn ' 



(4.16) 
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which is widely used in many low pressure steam studies is 
actually an incorrect deduction. For steam at low pressure % = 
1.3, Pr = 0.95, Gyarmathy's equation (with a = 0) should be 


a 


a 

c 


1 + 3 . 9 3Kn ' 


(4.17) 


The difference between the current model with a = 2 and 
Gyarmathy's model is shewn in Figure 4.3. It shows that while 
negligible error is incurred for both high and low values of Kn, 
significant deviation exists for moderate Kn . 


4.2 Droplet Temperature 

The droplet growth rate can also be found by molecular kinetic 
considerations [1.09], [1.19] as 
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With the modification which accounts for the vapor moving toward 
zhe droplet with a finite bulk velocity and the condensation 
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coefficient different from the evaporation coefficient [4.07], 
[4.03], the equation is presented by Young [1.05] as 
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By equating the growth rate found from kinetic theory Equation 
(4.18), or (4.19) with the general Equation (4.15), and 
eliminating dr/dt, we can obtain an equation for droplet 
temperature as done by Gyarmathy [1.09], Campbell and Bakhtar 
[4. OS], and Young (1.05]. However, the expressions are quite 
complicated, and the solution requires an iterative procedure. A 
much simpler equation proposed by Gyarmathy [1.09] can also yield 
the droplet temperature without the iteration procedure 
previously mentioned 


r 

T = T (p) - AT — (4.20) 

f s * r 

Young [1.05] argued that a net condensation is occurring and the 
droplet may not be in equilibrium with the vapor, yet Gyarmathy 's 
equation will only be valid by assuming the pressure of the 
droplet (after corrected for surface curvature) be the same as 


the vapor pressure. Nevertheless, it is shown by Gyarmathy [1.09] 
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that Equation (4.20) yields essentially the same result as the 

previous procedure discussed in the beginning of this paragraph. 

In addition , Equation (4.20) will reflect the physical meaning of 

the critical radius. If the droplet radius is grater than 

critical radius r > r , then T > T # the vapor will continue 

cr f g 

condensing on the droplet; If the droplet radius is less than 

critical radius r < r , then T < T , the droplet will 

cr f g 

evaporate. When r = r , T = T , the droplet size will remain 

cr f g 

the same as long as the vapor condition does not change. Also, 
the value of condensation coefficient and the relation between 
the condensation coefficient and evaporation coefficient remains 
unclear. Thus even with the complex iteration procedure, the 
final solution is still uncertain. Hence, Equation (4.20) will be 
used in the current model. 


4.3 Droplet Radius and Spectrum 

With Equation (4.2), (4.15) and (4.20), we have the complete 
information to calculate the droplet radius. Since droplets are 
continuously generated in different sizes by spontaneous 
nucleation during the expansion, we would divide the expansion 
process into many intervals, and the droplets generated in each 
interval will be treated as one group. Different groups of 



droplet then would have different numbers, sizes and 

temperatures, i.e. droplet spectrum. 3efore we can derive the 
equations for accounting of different groups, two questions have 
to be answered. First, will the droplets collide with each other 
thus affects the spectrum? Second, what is the effect of the 
foreign particles? 

4.3.1 coagulation 

Coagulation means the formation of larger droplet by collision of 

small droplets. Thermal (Brownian) motion is the main cause of 

the coagulation. Stein [4.10] has shown that for a typical 

2 

transit time through a condensation zone of 10 jisec, a droplet 

o 

of radius 30 A will undergo approximately one collision with 
another droplet, which is a frequency considered to be negligibly 
low to result in coagulation. Wu [4.11] also estimated and 
concluded that time for coagulation process is longer than the 
nucleation process, thus the effect of coagulation is negligible. 

4.3.2 effect of foreign particles 

Oswatitsch [4.12] first showed that the condensation by 
heterogeneous nucleation does not contribute significantly to the 


amount of condensate during rapid expansion in supersonic nozzle. 
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It is found that even by assigning generous growth laws on 
foreign nuclei, the time available is too short for the formation 
of appreciable amount of condensate. The number of nuclei 
generated by spontaneous nucleation far exceeds the number of 
existing foreign particles. In an interesting experiment, Ducklo 
and Pouring [4.13] artificially introduced a large number of 

3 

smoke particles (up to 10 particles) at the inlet of a 
supersonic nozzle. The result showed that there was no effect of 
seeding on the condensation process in the nozzle. It is worth 
mentioning here that condensation on the nozzle wall surface was 
estimated by Gyarmathy [4.04] to be negligible. 

4.3.3 equations for droplet accounting 

The radius of a group of droplets formed at location x and 

i 

currently located at location x can be found by integrating 
Equation (4.2) from x. to x 

a (x >; x') [T (x.jx 7 )-: (x 7 )] dx 7 
i i f , i i g 

pTh (T (x 7 ) ) -h (T “(xTx'Y) j“C(x'y 

f g g f f , i i 

(4.21) 

From Equation (4.15), heat transfer coefficient car. be expressed 
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a s 


a (x ; x) 
i i 


k / r (x ; x) 

g ' 11 

1 ~ 16 y ~Kn 

(1 + 4 Kn) + 3 y+lPr 


( 4 . 22 ) 


From Equation (4.20), the droplet temperature can be expreruiod a s 


T 

f, i 


(x. ; 
i 


x) 


T (x) + AT (x) 
9 


1- 


r {x, ) 

cr l 

r (x < ; x) 


(4.23) 


If we denote wetness as 


w . ( x . ; x ) 

i i 


4 3 

J (x . ) T tt r , (x. ; x) p_ 


i 

C (x. ) 


(4.24) 


, the average droplet temperature can be defined as 


T (x) 



w (x. ; x) 
1 1 


T . (X. ; x) dx. 
f , 1 1 1 


w (x) 


(4.25) 


, and the total wetness will be 



w. (x. ; x) 
i i 


dx 

i 


w(x) 


(4.26) 



CHAPTER 5 . 


MATHEMATICAL FORMULATION OF MODEL AND 


NUMERICAL SOLUTION TECHNIQUE 


A model of 

supersonic 

flow of a vapor 

in a Laval 

nozzle 

wa s 

developed . 

Coagulation 

, as discussed 

in chapter 

4 is 

not 

significant 

for nozzle 

flow and therefore 

was not considered, 

. It 


was also assumed that the effect of foreign nuclei is negligible 
in the flow. The model consists of the constitutive equations, 
r.ucleation theory equations and equations for droplet growth and 
account ing . 

5.1 Nucleation Rate 

The equations of nucleation rate were derived in chapter 3 and 
will not be presented here again. Equation (3.1) is for classical 
theory, and Equations (3.2), (3.3) and (3.4) are for Deich's 
theory. In the cases where expansion starts from superheated 
vapor, Equations (3.5) and (3.6) are used for Deich's theory. 
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5.2 Droplet Growth and Accounting 

For obtaining the droplet growth and accounting the number of 

droplets of various sizes, the equations for critical radius, 

r , (minimum radius of a droplet which can survive); drooled 

cr 

temperature, T ; droplet radius, r ; average droplet 

f , i i 

temperature, T ; heat transfer coefficient, « ; wetness, w ; 

f i i 

total wetness, w, are needed. These equations have been derived 

and discussed in chapter 2 and chapter 4. A summary of the 

equations to be used is presented as follows. 

Critical Radius r (x) : Equation (2.7) 

c r 

Droplet Radius r (x ; x) : Equation (4.21) 
i i 

Heat Transfer Coefficient a (x ; x) : Equation (4.22) 

i i 

Droplet Temperature T (x ; x) : Equation (4.23) 

Wetness w . ( x ; x) : Equation (4.24) 

Average Droplet Temperature T (x) : Equation (4.25) 

Total Wetness w(x) : Equation (4.26) 


5.3 Constitutive Equations 


Instead of using the differential form as in most investigations, 
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the integral form equations are used whenever possible. This 
greatly improves the computational efficiency. The equation for 
axial velocity, c, was obtained using the differential form of 
the continuity, momentum and energy equations. Detailed 

derivation can be found in Appendix D. It is assumed that the 
droplet travels with the same velocity as the vapor phase, and 
the vapor phase can be treated as an ideal gas. There are seme 
instances where real gas effect play an important role in the 
expansion process such as condensation of carbon-dioxide [1.19], 
high pressure steam [1.05], [5.01] etc. However, this effect is 

not important for low pressure steam as in the cases for 

experimental data used for comparison with the theox*y and the 
real gas effect can be included into the analysis and computation 
without difficulty whenever necessary. 

The temperature dependence of the various properties are 

considered in the analysis. The temperature range of interest for 

0 

the study cases in this work is between 0 to 130 C. The 

temperature dependent properties needed in the analysis include 

c , c , p , u, cr, and k . Among them, the constant value c = 

p f * 3 f 

4193 J/ (Kg • m) is used. Compared with c^ at the saturation line 
[5.01], this will give the maximum error < +0.5% between 0 — 100 
°C, and < -1.5% between 100 — 130 °C. The choice cf such value 
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is justifiable, for the higher temperature range is rarely 
encountered in the current cases studied. Also, since the actual 
c _ of the compressed liquid droplet tends to be lov;er than that 

o 

at the saturation line, the average error between 100 — 130 C 

will be even smaller. For similar reason and simplicity, the 

constant c = 1993 J/(Kg*m) is also used. Nevertheless, the 
P 

temperature dependence can be incorporated into the analysis 

without difficulty when different . temperature range is of 

interest. Surface tension cr(T ) is described by Equation (3.9) to 

g 

o 

an accuracy of ±0.2% between 0 — 13 0 C (compared with data 

adopted from [5.02]). Thermal conductivity k (T ) can be 

g g 

calculated to an accuracy of ±0.5% in the temperature range 0 — 
130 °C (compared with data adopted from [5.03]) as 

k (T ) = 1.82 x 10' 2 + 5.8343014 x 10~ 5 T 

g g g 

- 5.888523 x 10‘ 8 T 2 - 4.51 x 10" U T 3 (5.1) 

g g K-m 

where T is in °C. Viscosity u(T) can be calculated to an 

g g 

0 

accuracy of ±0.4% in the temperature range 0 — 130 C (compared 
with data adopted from [5.03]) as 

u(T ) = (0.393 T + 30.2) x 10 

g g 


Pa • s 


(5.2) 
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o 

where T is in C. The liquid phase density p (T ) can be 
g t t 

expressed to an accuracy of ±0.1% in the temperature range 0 — 
130 C (compared with data adopted from [5.03]) by the following 
equat ion [5.04] . 


P,( V = P c 1 + E D . ( 1 - =r ) 


1=1 


"1 


where 


p = 317.0 Kg/m 3 

T = 647.236 K 
c 

D = 3.6711257 

1 

D = -2 .8512396 x 10 1 
2 

D = 2.2265240 x 10 2 

3 

D = -8.8243352 x 10 2 

4 

D = 2.0002765 x 10 3 

5 

D = -2.6122557 x 10 3 

6 

D = 1.8297674 x 10 3 

7 

D = -5.3350520 x 10 2 

8 


(5.3) 


The constitutive equations are listed for both single component 
condensing vapor and two component (condensing vapor and inert 


gas) models. 
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5.3.1 single component model 


Continuity p : 

g 


P 


g 


c A 
in in 

c A 


p, (1-w) 
in 


( 5 . 4 ) 


Equation of velocity c : 


dc 

dx 



dA 1 
dx A 


wc dT c T -h' , 

t t p , g g t , g dw 

(l-w)c T dx (l-w)c T dx 
p, g g p , g g 


( 5 . 5 ) 


where 


M 


2 


(c - R)c 2 
P 


c RT 
p g 


h' = h (T ) 
f ,g g g 


h (T ) 
f f 


The derivative of total wetness 


dw 

dx 


in Equation 


( 5 . 5 ) 


is derived 


in Appendix C to be 



4 5 


dw 14 3 . 

-T- X) = — — r T[ r (X) p f J(X) 

dx c(x) 3 cr f 




4tt 


c (x) 


y. =0 

i 


r ( x ; x) a {x ; x) J ( x . ) [ T . ( x . ; x) - T { x ) ] dx. 
ii ii i f , l l g i 

[h (T (x)) - h (T . (x.; x))] c(x.) 

9 3 f f , i : i 


(5.6) 


Energy equation T : 
9 


2 

%- + (l-w) h (T ) + Wh (T ) = h 
2 g g f f 0 


(5.7) 


Equation of state p, (for vapor and supercooled region) : 


P = P n R T a 

9 9 


(5.8) 


5.3.2 two component (condensing vapor and noncondensing gas) 
model 

The two components considered here are the condensing vapor and 

the noncondensible gas which does not react with the condensing 

vapor, nor has condensation. It is assumed that vapor phase 

temperature for both components is the same and is represented as 

T . The pressure of each component is represented by partial 

g 



4 6 

pressure. The mass fraction of the condensing vapor contained in 
the mixture at the beginning of the expansion is denoted by . 
If subscripts n and v denote the noncondensible gas and the 
condensing vapor respectively, partial pressure can be defined as 


P = P + P 
v n 


( 5 . 9 ) 


p = y p 

V 


( 5 . 10 ) 


p_ = (1 - y) p 


( 5 . 11 ) 


(w - w) / M 
0 v 

(1 - w ) / M + (w - w)/ M 
0 ' n 0 1 v 


( 5 . 12 ) 


For two component model, the nucleation rate and droplet growth 
and accounting equations are the same as single component model 
except that partial pressure p has to be used, and all the 
variables should now be referring to the condensing vapor in the 
two component case e.g. gas constant R, number of molecules per 
unit mass N etc. The supercooling and supersaturation will now be 
evaluated as 


AT = T (p ) - T 
s v g 


( 5 . 13 ) 
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A = In (p / p (T )) (5.14) 

' v ' s g ' 

The expressions of the constitutive equations are, however, 
somewhat different from the single component model. They are 
presented as follows. 

Continuity p : 




A 


1 - w 


CA 


(5.15) 


Here p is the density of the whole vapor phase including 

g 

condensing vapor and inert gas. A mean molecular weight can be 
defined as (see Appendix 2 for derivation) 


M 

g 


1 - w 

(1 - w ) / ~M + (w - w) / M 
0 1 n 0 7 v 


(5.16) 


The relationship between p and p , p can be shown to be 

g v n 


(Appendix E) 
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M 

g 



( 5 . 17 ) 


Equation of velocity c : 

The equation of velocity can be derived with the same procedure 
as one component case (explained in Appendix D) 


dc c I 1 dA 

** = m * 2 - i r ^ 


1 / M 

' V 

(1 - w )/ M + " ( w - w)/ M 
0 7 n 0 7 v 


( 5 . 18 ) 


(c - wc ) T 
pO pv g- 


dw 
dx + 


wc dT \ 

£ il 

(c - wc ) T dx 
pO pv g J 


where 


«* 2 - r f 

a v 


RU1 - w )/ M + {w - w) / M ] 
0 7 n 0 7 v 


c - wc 

pO pv 


h 


f V 


= h 

V 


(T ) - h (T ) 
g £ f 


c 

pO 


( 1 - w ) c + w c 

0 pn 0 pv 
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Energy equation T 

g 


C 

2 ~ 


+ c (T 
po q 


T ) - w h 
0 fv 


( 5 . 19 ) 


Equation of state p : 


P 



w ) / M +(w - w) / M 

0 7 n 0 7 v 

1 - w 


T 

g 


(5.20) 


5.4 Numerical Solution Technique 

The solution procedure used for solving these equations was as 

follows. Initial values of the temperatures, velocity c, density, 

etc., were knowmi . Finite difference approximations were obtained 

1 dA 

for d(AT)/dx and — — . The solution procedure consisted of 

A dx 

guessing the values of T , pressure p, velocity c, and radius r , 

g i 

and using an iterative procedure to solve all the equations to 
get a new set of the values of the variables at the grid point. 
The procedure was then repeated for the next point . A fourth 
order Runge Kutta method was used for solving Equations (5.5) and 


(5.18) . 
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It can be observed that Equations (5.5) or (5.18) will break dow n 
v.'hen the value of the term M/(l - w) or (M - 1) equals to .1 . 
Numerically, problem arises even when M~/{ 1 - w) or (M - 1) is 
very close to 1. It is observed from Equations (5.5) and ( 5 . 18 ; 
that if wetness is not zero, this singular point would occur at a 
distance upstream of the actual nozzle throat . Hew large this 
distance is would depend on how large the wetness is. The 
r.ucleaticn usually starts somewhere upstream of the nozzle 
throat, hence, c he wetness is usually greater than zero. However, 
the nucleation rate is usually not high there, and tne wetness is 
usually small, consequent ly , this distance is usually small as 
well . 

The singularity characteristic in two-phase flow equations was 
discussed by many investigators [5.05], [5.06], [5.07]. Due to 
this singular nature, a pure numerical scheme, if employed for 
solution, will either yield erroneous results or will become 
unstable and fail. Therefore, a remedial procedure must be 
adopted in this region. One way is to use the differential form 
zi the momentum equation. The differential form of the momentum 
equation can be expressed as 
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dc 1 - w dp 

c _ — ^ -r— =0 , . 

ax p ax 

5 

With proper finite difference approximation (Appendix D) , the 
velocity at next grid point 2 can be computed from the corn' 1 * 
grid point 1 as 


c 

2 



(w + w ) - 4 
1 2 


+ P 2 


(Pn 


P. ) 

i 


(5.22) 


Although the above mentioned method is easier to implement than 
using Runge-Kutta method on the equation of velocity, its 
accuracy is not as good. Fortunately, the axial range of nozzle 
where M 2 / (1 - w) or (M +2 - 1) is very close to 1 is very short. 
Also, this usually occurs very' close to the nozzle throat where 
nucleation rate and wetness are small. Therefore, this remedial 
procedure does not affect the final result. 

In some of the cases studied, the aforementioned technique either 
will not converge or dees not give the physically acceptable 
supersonic branch of solution. Under this circumstance, another 
technique has to be used. As we noticed that the pressure ratio 
always starts to deviate from the isentropic pressure ratio at 
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quite a distance downstream of the the singular point. Therefore, 
the pressure ratio can be taken as the isentropic pressure ratio 
at the singular point. The velocity is calculated utilizing t ho 
differential form of the momentum equation. Until the "dangerous 
zone" is passed, the computation is switched back to equation oil 
velocity. Depending on the different cases studied, the length of 
this zone varies. Different machine accuracies associated with 
the different computers used do not make appreciable change on 
the length of this zone. Instead of using the isentropic 
pressure, the extrapolation of the properties can also be adopted 
to solve the singularity problem. When doing so, one has to make 
certain that the solution after extrapolation falls into the 
correct quadrant of the solution space. Usually, a criteria 5p < 
0 and <5x > 0 is used [5.05] where Sx is the distance for 
extrapolation. In addition, assumptions of negligible nucleation 
rate and negligible change of droplet temperature between the 
extrapolation points have to be made, and a simplified droplet 
growth rate equation has to be applied. Some other techniques are 


discussed 

by Gyarmathy [1.09]. 

Some 

studies neglect 

the 

nucleation 

effect entirely before 

the 

throat and start 

the 


computation from the throat, thus avoiding the singularity point 


[5.03] . 



CHAPTER 6 RESULTS AND DISCUSSION 


6.1 Single Component Results 


6.1.1 Pressure ratio and mean droplet radius 


Although the two nucleation theories are substantially different, 


both 

need the value 

c f the 

surface tension 

for 

droplets with 

smal 1 

radius. The 

surface 

tension values 

for 

small radius 

droclets are net available. 

According to 

many 

investigators 


[3.09, 3.10, 3.11, 3.12, 6.01] the surface tension of a small 

droplet is different from that of a flat surface and one is net 

sure whether it decreases or increases. Therefore, various 

values of the surface tension ratio (c r/<r ) were tried to obtain a 

00 

good agreement between the calculated pressure ratio with the 
experimental results. Sixteen test cases for which experimental 
data were available were studied. A summary of inlet total 
conditions, Wilson pressure, surface tension ratio, the mean 
radius as well as the standard deviation at the end of 
computation (not necessarily the nozzle exit, depending upon the 
nozzle geometry information available) is presented in Table 6.1. 
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The pressure ratio along the nozzle length obtained using 
classical and Deich's nucleation theories compared reasonably 
well with the experimental data. Results of pressure ratio and 
droplet size distribution for all the sixteen cases are provided 
in Figures 6.1 to 6.43. It is observed that in all cases, the 
final droplet size distribution is quite monodisperse . However, 
Deich's theory gives slightly better . agreement with experiments 
in general. Computation also shows that Deich's theory gives 
better monodisper s ion regarding the droplet size distribution and 
smaller mean radius (see Figure 6.6). 

Experimental information on mean droplet size are scarce and only 
the data provided by Gyarmathy and Meyer [1.13] was available for 
comparison. Figure 6.49 shows the comparison between the mean 
radius predicted by using Deich and classical nucleation theories 
and the experiment (Gyarmathy and Meyer Nozzle data, Table 6.1). 
The mean droplet size obtained by using Deich's theory is 
reasonably close to the experimental data, while the classical 
theory tends to overestimate the droplet size. However, it 
should be noted that as mentioned earlier the surface tension 
ratio were varied to obtain good agreement with the experimental 
pressure ratio and the same value of the surface tension ratio 
was used for ascertaining the droplet sizes. 
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6.1.2 Typical results of some important variables 

The results of the sixteen test cases investigated showed si.no Jar 
trends. Therefore, we will only present the result of 

Barschdorff nozzle test no. 7 as representative of the data. 
Figure 6.50 show's the variations - of nucleaticn flux and 
supercooling along the nozzle axial distance. Nucleation starts 
after a certain ajnount cf supercooling (in this case about 17 C) 
is reached. As steam expands, nucleation reaches its peak 
slightly upstream of the peak cf supercooling and drops to zero 
suddenly after the supercooling peak where spontaneous 
condensation occurs due to large number of nuclei. Since most of 
the nucleation occurs around the nucleation flux peak, this 
accounts for the monodisperse nature of droplet size. The latent 
heat released due to this spontaneous condensation causes a sharp 
pressure increase, as showu in Figure 6.1 through 6.48. Note that 
pressure rise occurs at the peak of supercooling cr right after 
the peak of nucleation rate, which is the start of the 
spontaneous condensation and is defined as the Wilson point. The 
latent heat also causes a sudden increase of the vapor 
temperature (Figure 6.51) and entropy increase (Figure 6.52), and 
results in bringing the vapor toward an equilibrium state rapidly 
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from a supercooled state. Consequently, supercooling decree s -v- 
rapidly during spontaneous condensation. At the end o>; 

spontaneous condensation, usually defined at where supercooling 
stops to decrease and remain constant, the state of vapor is 
close enough to the equilibrium state (about 3 C supercooling >■.= 
this typical example), and the expansion proceeds almost like 
normal equilibrium expansion. The pressure (Figure 6.1 through 
6.43) and vapor temperature (Figure 6.51) decrease again similar 
to where normal equilibrium expansion takes place, and the 
entropy increases in a much slower pace (Figure 6.52). 

Figure 6.51 shows the vapor temperature T , liquid temperature 
T and the eoui librium saturation temperature at the vapor 

f 

pressure T (p) , along the nozzle axis. The result indicates that 
s 

the liquid temoerature T is the same as T at the beginning of 
nucleation. At the end of spontaneous condensation, T and T^ 
are very close, which means that the liquid phase is almost in 
equilibrium, while the vapor phase remains at a constant 

supercooling (about 3 C) . We can also observe from this figure 
that wetness becomes appreciable after the sudden pressure rise. 
Before the Wilson point is reached (i.e., sudden pressure rise or 
start of spontaneous condensation), the increase in the wetness 
is mainly due to the new nucleation which is so small that the 
steam is practically dry. After the Wilson point, which is also 
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the end of nucleation and the start of spontaneous condensation, 
the existing large amount of nuclei (nucleated around, the 
nucleation flux peak) provide a large surface area for 
condensation which causes a rapid increase in wetness. 

Mach number and entropy are shown in Figure 6.52. The variation, 
in the Mach number and the entropy along the nozzle axis is 
similar to the pressure variation. • Before the occurrence of 
spontaneous condensation, Mach number increases, entropy remains 
constant and expansion is i sent rep ic . During spontaneous 
condensation, Mach number decreases and entropy increases because 
of the pressure rise and release of latent heat. After the end 
of spontaneous condensation the vapor is almost in a 
thermodynamic equilibrium state, the expansion is nearly 
isentropic and the Mach number increases again. 

Figure 6.53 shows the growth of mean radius and Sauter mean 
radius of the droplets and the variation in critical radius along 
the nozzle axis. Mean radius is defined with the same amount of 
wetness and the same total number of droplets. Sauter mean radius 
is defined by the ratio of the total volume to the total area of 
all the droplets. Sauter mean is slightly larger than the mean 
radius in all the study cases here. It should be noted that the 
smaller the difference between the two mean, the better the 
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monodispersion would be. The droplet growth rate is basically 
proportional to the temperature difference between the droplet 
and supercooled vapor. In the nucleation zone (before the advent 
of spontaneous condensation) in spite of the increasing 
temperature difference (Figure 6.51), the droplet growth xi>' <• 
remains small due to the presence of the large number of newJy 
formed small nuclei. At the end of nucleazion and at the 
beginning of the spontaneous condensation, mean radius grows 
rapidly because no new small nuclei are being formed. Near the 
end of spontaneous condensation, the temperature difference 
becomes small, therefore, the growth rate slews down again. The 
critical radius decreases in the nucleation zone, which also 
justifies the consideration of the second component of nucleation 
in Deich's nucleation theory. The results of droplet size 
distribution at the end of the expansion are presented in Figure 
6.1 to 6.48. From the figures, we can see that good 
monodispersion exists for all study cases. 

The nozzle expansion rate and velocity for Barschdorff No. 7 are 
shown in Figure 6.54. The shape of the expansion rate curve of 
each case studied mainly depends on the nozzle geometry. The 
velocity curve displays similar trend as the Mach Number. 


Figure 6.55 depicts the Knudsen Number and Prandtl Number. In all 
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the cases studied, Prandtl Number decreases before Wilson point 
is reached, rapidly increases during the spontaneous 
condensation, and stabilizes to remain nearly constant after 1 -ho 
spontaneous condensation. The change in Prandtl Number is small, 
and the value is around .95 for all study cases. The 
Number always increases at the early stage of the expansion where 
the pressure is decreasing and the droplet growth rate is very 
slow. When droplet growth rate starts to pick up, the value of 
Knudsen Number drops sharply until the droplet growth rate slows 
down. In most cases the effects of pressure drop and droplet 
growth eventually balance cut, and the Knudsen Number would 
remain fairly constant through the rest of the expansion. The 
range of the value of Knudsen number varies nozzle by nozzle. In 
the cases studied, it ranges from one to a few hundreds. Thus 
we know that for most cases the flow is under free molecular flow 
regime (Kn > 3) and sometimes under slip and transition flow 
regime (0.01 < Kn < 3). The continuous flow regime (Kn < 0.01) 
is almost nonexistent. 

The two components of the nucleation rate for Deich's nucleation 
theory are presented on Figure 6.56. The second component has 
always been a few order of magnitude smaller than the first 
component in all our study cases. However, depending on the 
condensing vapor used, expansion rate, and the initial condensing 
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vapor content, it is possible that the second component may 
overwhelm the first component. 


6.1.3 Correlation for surface tension ratio 

In this study the surface tension ratio values were varied to 
obtain agreement between the theoretical and experimental 
pressure ratios. It is therefore important to establish relevant 
correlation fcr the surface tension ratio so that a proper value 
can be used fcr designing new nozzles and analyzing nozzles fcr 
which experimental data does not exist. Thus, obtaining 
correlation for surface tension ratio for the nucleation theories 
is very important. Various correlations were tried and the 
results showing inlet total temperature, Wilson pressure, inlet 
total pressure and inlet total entropy for the two nucleation 
theories are presented in Figure 6.57 to 6.64. Apparently, inlet 
total temperature does net yield acceptable correlation for both 
nucleation theories. The correlation yielded by Wilson pressure 
seem to show a positive trend. However, since the Wilson pressure 
is not known a pnJuonL, a correct value of the surface tension 
ratio cannot be assigned. A more difficult iterative procedure is 
needed to implement the Wilson pressure correlation in nozzle 
design. On the other hand, the usage of inlet total entropy and 



pressure correlations are rather straightforward and need no 
iteration procedure since the inlet conditions are always known . 
In addition, the Wilson pressure correlation data presents the 
larger scatter than inlet total pressure and inlet total entropy 
correlations . 

The correlations utilizing inlet total pressure and inlet total 
entropy seem to be the most promising. However, the slope of the 
correlation based upon inlet total pressure at the region of 
small total pressure is so steep that a small uncertainty in the 
correlation would trigger a large error in the predicted surface 
tension ratio. As a result, the correlation between the inlet 
total entropy and surface tension ratio is chosen. From Figure 
6.63 and 6.64 we found that the trends on the variation of 
surface tension ratio with inlet total entropy for Deich's and 
classical nucleation theories are similar. However the scatter 
for the classical theory data is larger than that for Deich's 
theory. The correlation coefficient (which is defined in Appendix 
G) for Deich's theory can be computed as -0.9394 and for 
classical theory to be -0.9337. The larger correlation 
coefficient for Deich's theory shows that it gives better 
correlation. In the next section it will be shown that the 
pressure ratio, mean droplet size and droplet size distribution 
are quite sensitive to the value of the surface tension ratio. 
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Therefore, the correlation for Deich's theory with less scatter 
was adopted. Additionally, the classical theory, as indicated 
earlier, tends to overestimate the droplet size. 


6.1.4 Effects cf parameters 
A. effect of surface tension ratio 

The effect of the variation of surface tension ratio on pressure 
rise associated with spontaneous condensation, mean radius and 
droplet size distributions are shown in Figures 6.65, 6.66 and 
6.67 respectively. These results were obtained by using different 
values of surface tension ratio for the Barschdorff nozzle No. 7 
for the same inlet conditions p^ = 78323 Pa, T q = 388.44 K . 
Higher surface tension ratio which gives a smaller nucleation 
rate and broader nucleation zone will delay the onset of 
spontaneous condensation {indicated by pressure rise) thus 
decreasing the Wilson pressure. The total amount of heat transfer 
can be roughly measured by the extent of the pressure rise. From 
Figure 6.65 we found that the total amount of heat transfer is 
about the same for different surface tension ratios. The mean 
droplet radius is smaller at first in the case of larger surface 
tension ratio since the nucleation starts later compared with the 
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case of smaller surface tension ratio. When spontaneous 
condensation takes place later, the mean radius starts to grow 
rapidly. With the approximately same amount of total heat 
transfer and the smaller number of nuclei given by larger surface 
tension ratio, more heat transfer takes place on each individual 
droplet. Consequently, the mean droplet radius of larger surface 
tension ratio eventually outgrows those of smaller surface 
tension ratio (see Figure 6.66), and causes a slightly larger 
variation in droplet radius (Figure 6.67). The results are 
relatively very sensitive to the change in surface tension ratio. 
This reaffirms our previous statement concerning the inability to 
obtain an acceptable correlation for surface tension ratio using 
classical nucleation theory due to relatively large scatter in 
correlation data. 


E. effect of inlet total temperature 

The effect of inlet total temperature is investigated by 
comparing Barschdorf f ' s nozzle test No. 4, 7 and 10 at an inlet 
total pressure of 11.36 psia. From Figures 6.63, 6.69 and 6.70, 
it is observed that for a given nozzle geometry and inlet total 
pressure, a higher value of inlet total temperature delays the 


onset of spontaneous condensation and therefore decreases the 



Wilson pressure and mean droplet size and gives a relatively 
smaller variation in droplet radius. The inlet total temperature 
in general does not affect the nucleation rate. Thus, the mean 
droplet radius for higher inlet total temperature does not have 
the outgrown effect similar to the cases of different surface 
tension ratio. 


C . ef fect of inlet total pressure 

The effect of inlet total pressure is obtained by coirparing the 
results for the tests performed on nozzle of Barschdorff with 
total temperature T q = 233. 44K, and different pressure at 75974 
Pa, 78323 Pa and 30673 Pa. Proper surface tension ratios are used 
according to the correlation developed in the previous section. 
The inlet total pressure like the inlet total temperature in 
general does not affect the nucleation rate. It is observed from 
Figures 6.71, 6.72 and 6.73 that lower inlet total pressure 
delays the onset of spontaneous condensation and therefore 
decreases Wilson pressure and mean droplet size and provides a 
relatively smaller variation in the droplet radius. The effect is 
similar to that of higher inlet total temperature. However, the 
results are more sensitive to the change of inlet total pressure 


than to the change of inlet total temperature. 
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D. effect of expansion rate 


In the previous sections, the effects of parameters are 
considered for the same nozzle with different inlet total 
conditions and/or different surface tension ratio. The expansion 
rate is determined by the nozzle geometry. Hence, in this 
section we shall consider the effects of different expansion 
rates (thus different nozzle geometries) with the same inlet 
total conditions and surface tension ratio. The results obtained 
for the three nozzles designed in the next chapter are used i.e. 
Nozzle NZ11 with expansion rate 2500 1/sec, NZ12 with expansion 
rate 5000 1/sec and NZ13 with expansion rate 10000 1/sec. The 
inlet condition, total temperature 500 K and total pressure 1.7 
MPa is used. 

In Figure 6.74 we can see that for higher expansion rate the 
larger area ratio causes smaller pressure ratio and the onset of 
condensation takes place earlier. The higher expansion rate 
generates higher nucleation rate and narrower nucleation zone. 
Although the total amount of pressure rise, thus the total amount 
of heat transfer, is about the same for different expansion rate, 



the pressure rise is steeper for the higher expansion rate case 
since heat transfer takes place faster due to more nuclei 
available. The mean radius of lower expansion rate nozzle is 
smaller compared with the higher expansion rate nozzle due to 
delay of the onset of condensation. However, when spontaneous 
condensation occurs later, with the same amount of total heat 
transfer and less droplets, more heat transfer takes place for 
each individual droplet in the lower . expans ion rate case. As a 
result, the mean radius of lower expansion rate nozzle eventually 
outgrows the one of higher expansion rate nozzle. 

Figure 6.75 shows the droplet sizes at the end of the expansion 
through the nozzle and we can observe that the higher expansion 
rate nozzle NZ13 yields smaller mean radius =* .1 fim as compared 
to .25 pm for the nozzle of an expansion rate 2500 sec . The 
variation in the droplet radius is also smaller for the higher 
expansion rate (Figure 6.76). 

From the preceding discussions, we can conclude that the droplet- 
size distribution and the droplet mean radius is a function of 
the inlet total pressure and temperature and the nozzle expansion 
rate. It is also sensitive to the value selected of the surface 
tension ratio. However, the variations (fluctuations) in the 
droplet radius are relatively insensitive to these parameters. 
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The variations in the droplet radius, for most of the cases 
studied, were within ± 0.01 pm. This indicates that the droplet 


sizes stay 

relatively monodisperse 

for various 

values 

of 

the 

parameters . 

A 

summary of 

the mean radius 

at 

the 

end 

of 

ccmputat ion 

for 

the sixteen 

cases is 

provided 

in 

Table 

6.1. 

The 

droplet radius 

varies from 

0.003 pm 

to 0.074 

pm 

depending 

upon 


the inlet total condition and the expansion rate. We can thus see 
the importance of manipulating the parameters to design a nozzle 
which can provide monodisperse droplets of desired radius. 


6.2 Two Component Results 


There are very few experimental data available for two component 
(one condensing vapor and one noncondensing vapor) spontaneous 
condensation. The vapor temperatures of the two available 
experiments both go well below the sublimation line [6.02], 
[6.03]. The uncertainties on the free surface energy, the 
spontaneous crystallization rate, the thermodynamic properties of 
supercooled liquid and ice nuclei, and the validity of the 
nucleation theories etc. makes the comparison between the 
theoretical numerical result and experimental data hardly 
meaningful. However, in order to investigate the effect of 
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initial vapor content on the spontaneous condensation, the nozzle 
of Barschdorff No. 7 is used in the numerical computations. 


6.2.1 Effect of some important variables cn the results 

According to gasdynamics, we know that isentropic pressure ratio 
is not only function of area ratio (i.e. nozzle geometry), but 
also function of ratio of specific heat y. For a two component 
vapor, the ratio of specific heat y depends on the initial vapor 
content. As a result, the pressure ratios for different initial 
vapor content for the same inlet condition shown in Figure 6.77 
are different. Some other typical results for different initial 
vapor contents are shown in Figure 6.73 to 6.54. Comparing the 
results for wO=l with single component result shown in Figure 
6.50 to 6.56, we find that the results are identical. This can be 
seen as a validation of the theory and the numerical program. 


6.2.2 Effect of initial vapor content 

The effect of initial vapor content has a basic difference from 
those parameters presented in the single component case. From 
Figure 6.77 we notice that smaller initial vapor content gives a 
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smaller pressure rise. This is because smaller heat transfer due 
to phase change occurs due to less condensing vapor available. In 
single component case, all the parameters studied basically has 
no effect on the total amount of heat transfer. Smaller initial 
vapor content also delays the onset cf the spontaneous 
condensation indicated by pressure rise and in turn the Wilson 
pressure. However, contrary to higher surface tension ratio and 
lower expansion rate in single component case, smaller initial 
vapor content gives higher nucleation rate and narrower 
nucle&ticn zone. With smaller total amount cf heat transfer and 
more nuclei, less heat transfer takes place on each individual 
droplet. Therefore, smaller mean radius (Figure 6.73), better 
monodispersion cf droplet size (Figure 6.79), and less variation 
on supercooling (Figure 6.30 to 6.82) were found in smaller 
initial vapor content case. This in turn makes the mean radius 
growth rate slow (Figure 6.73) and generates less condensate and 
therefore, results in decreases in the wetness (Figure 6.89 to 
6.91) . 


From Figure 6.33 to 6.35 we can also observe that while Pr tends 
to be smaller for the smaller initial vapor content, it still 
remains in the range of 0.9 to 1.0. On the other hand, Kn changes 
substantially with changes in initial vapor contents. The smaller 
the initial vapor content the larger Kn would be. Consequently, 
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the flow tends to behave as a free molecular flow (Kn > 3). 
Another interesting phenomena found from Figure 6.86 to 6.83 is 
that the order of magnitude of the two components of nucleation 
flux of Deich's nucleation theory are closer in the case of small 
initial vapor content. This shows that the stable nuclei 
initially formed are so small in the low initial vapor content 
case that the unsteady effect (represented by the second 
component of Deich's theory) plays a stronger role. Figure 6.92 
to 6.94 shew that with decreasing initial vapor content, the 
actual velocity as well as the isentropic expansion rate 
decrease . 



CHAPTER 7 NOZZLE DESIGN FOR OBTAINING MONODISPERSE DROPLETS 


The objective of nozzle geometry design is to obtain proper 
A/A (x) such that the desired droplet mean radius can be 
generated for certain inlet total conditions. In the preceding 
chapter, we have discussed the effects of different parameters 

T d , and p on the mean radius of the droplets. We have also 

0 “ 0 

discussed the correlation bet ween <r/cr^ anc ^ ^he inlet total 
condition. These results form the basis of nozzle geometry 
design . 

7 . 1 Equations 

The relation ajnong T , p/p (x) , p(x), r( T(x)) f and c (T(x)) can 

0 0 P 

be expressed as [1.09] : 
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(7.1) 


in which P/P Q (x) 


is directly related to A/A (x) and y(T(x)) by 


A_ 

* 

A 



l 

1 




(y-u /y 


(7.2) 


Since T(x) is not known a priori, the coupling among the 

variables in the above equations is rather complicated. In most 

cases, however, y and c are constant in the temperature range of 

p 

interest. Thus for given T q and p(x), we can derive p/p Q (x) from 

* 

Equation (7.1), and then derive A/A (x) from Equation (7.2). A 
convenient way which does not violate the generality is to 
further assume a constant expansion rate. With the above 
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assumption. Equation (7.1) can be solved for p/p^(x) easily using 
Newton's Method. 


7.2 Design Procedure 


The following procedure is based on constant y and c^, and 

utilizes constant expansion rate and as design parameters . 

This procedure will lead to a preliminary nozzle profile. 

Step 1 s Choose certain T^ and p, solve for p/p Q (x) from 

Equation (7.1) and then solve fer A/ A (x) from 

Equat ion (7.2) . 

Step 2 : By changing p ( P ^ P A1 . . ) # we can obtain 

‘ 1 o 0 o limit 

different droplet mean radius : 

p T > r T 

*0 1 m 1 

P 0 1 * r m 1 

p is the highest pressure allowable for 

ulimit 

maintaining as gaseous phase at the chosen T q in 

Step 1. p greater than p is not allowed. 

v 0 Olimit 

Two options are available in this step: 

(a) If r is too small, follow Step 3; 
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(b) If r is too large, follow Step 5. 

Step 3 : There are two options in this step: 

(a) Start all over from Step 1, use the 

nozzle with the same T and lower p to 

0 

♦ 

obtain p/p (x) and A/A (x) ; 

0 

(b) Start all over from Step 1, use the 
nozzle with the same p and higher to 
obtain p/p^(x) and A/A (x) . 

Although both options lead to larger r_, there are 
some differences between the two options. Using 
option (a) one does not have to increase p^ 
substantially higher to achieve the same increase 
in r compared to option (b) . However, low 
expansion rate nozzles are relatively more 
difficult to fabricate due to the sina ) ] t\i 
variation in area ratio. With option (b) , 

depending on the type of expanding vapor and the 
mean droplet size desired, one may have to raise 
p^ to an extremely high level which may be 
difficult to implement experimentally. 

Repeat Step 2 until the desired r is obtained. 


Step 4 : 
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Step 5 : The two options exactly opposite to Step 3 are as 

follows : 

(a) Start all over from Step 1, use the 

nozzle with the same and higher p to 

obtain p/p (x) and A/A (x) ; 

0 

(b) Start all over from Step 1, use the 

nozzle with the same p and lower to 

obtain p/p (x) and A/A (x) . 

o 

Step 6 : Repeat Step 2 until the desired r is obtained. 

Variable expansion rate can be easily implemented without 
technical difficulties by the same procedure if the resulting 
profile is not satisfactory. The experience, trial and error, and 
state of art selection of p(x) will then have to come into the 
play. Even in the cases where y and c^ vary substantially with 
temperature change, this preliminary nozzle profile can be used 
as the prototype for further modifications. An iterative design 
procedure can be adopted to perform the modification to 
satisfaction. First run the numerical test on the preliminary 
nozzle profile to obtain temperature T(x) . Incorporate this T(x) 

into Equation (7.1) and (7.2), we can then solve for p/p^(x) and 

* 

the new nozzle profile A/A (x) and proceed with the design steps 
described above. Thus iteratively we can modify the design to our 



satisfaction . 


7.3 Some Examples Of Nozzle Geometry Design 

With combinations of different total- temperature and expansion 
rate, ficjure 7.1, 7.2 and 7.3 show the area ratio of a few sample 
nozzles designed for steam. By comparing the figures, v;e can see 
that for the same expansion rate, higher total temperature would 
give flatter nozzle shape. The same is true for the same total 
temperature with lower expansion rate. It is noted that the 
result of these nozzles are used to show the effect of expansion 
rate on some important variables in the previous chapter . 

i 

Using different total pressure, the mean radius at the exit of 
each nozzle is calculated and summarized in table 7.1. Such a 
table can serve as a design guideline and would be very useful to 
generate the geometry of the nozzle that could give the desired 
mean droplet radius. For example, if a mean radius of 0.16 pm is 
wanted, by looking at the table we could immediately find that 
for such nozzle, falls between 500 K and 425 K, expansion rate 
falls between 2500 1/sec and 10000 1/sec. Or one may use nozzle 
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NZ12 and fine tune the inlet total pressure between 1.7 M?a and 
.37 MPa to obtain the desired mean radius. 

Although different nozzles may all be able to achieve the same 
goal with proper inlet total pressure and/or proper length, not 
all the implementations are feasible. Tv; o determining factors 
are the degree cf slanting of the physical nozzle shape and the 
exit pressure ratio cf the nozzle. Highly slanted nozzle would 
have strong boundary layer effect, multidimensional effect and 
flow separation problems. Flat shaped nozzle presents less 
boundary layer effect and less flow separation problem, however 
creates more difficulties in precision and fabricating the 
nozzle. High exit pressure ratio would require either extremely 
high inlet pressure or the extremely high exit vacuum which are 
in general difficult to establish, with such extremely high 
pressure ratio, the strength cf the nozzle material will a. Iso be 
a limiting factor. Thus a compromise between all the limiting 


factors has to be reached for a feasible design. 



CHAPTER 8 CONCLUSIONS 


An analytical model to study spontaneous condensation of single 
component and two component (condensible vapor and noncondensible 
gas mixture) was developed. Two nucleation theories one based 
upon classical theory (Volmer and Frenkel), another Deich's 
theory were used. The following conclusions result from this 
study . 

1. Computer codes to solve the models for both single component 
flow and two component (one condensing vapor and one 
noncondensing vapor) flow in a condensing nozzle were developed 
and validated using available experimental results. The 
comparison between experimental and numerical results was good. 
The code provides the mean droplet size and number density, 
pressure ratio and quality as a function c: nozzle geometry and 
nozzle inlet conditions. 

2. The pressure ratio obtained from both classical theory and 
Deich's theory compared favorably with the experimental data, 
however, the results of Deich's theory have slightly smaller 
standard deviation. In designing nozzles for obtaining relatively 
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monodisperse droplets, Deich's theory should be used based on the 
following reasons: (a) No acceptable correlation for surface 
tension ratio can be established for classical theoiy; (b) 
Classical theory tends to overestimate the mean droplet radius; 
and (c) Deich's theory is expected to perform better on high 
expansion rate nozzles than classical theory because it also 
accounts for the nozzle expansion rate, which plays an important 
role in the spontaneous condensat ion , in the nucleation process. 

3. Surface tension ratio can be correlated with inlet total 
entropy for use with Deich's theory. A correlation for steam flow 
was developed and used success ful ly . 

4. Lower inlet total pressure and higher inlet total temperature 
lead to a delay in pressure rise, smaller mean radius of the 
droplets after spontaneous condensation, and smaller variation .in 
droplet sizes. 

5. Higher expansion rate will generate higher nucleation rate 
and have steeper and earlier pressure rise. It also results in 
smaller mean radius and better monodispersion. 

6. The influence of all parameters on monodispersion is rather 
small indicating that the spontaneous condensation process will 



30 


always result in relatively monodisperse droplets. 

7. The influence of noncondensing vapor was investigated ujj.ng 
the two component flow model. Lew initial condensing vapor 
content results in a later and flatter pressure rise, smaller 
mean radius, larger number of droplets and better mcnodisper s ion . 
It also results in higher Kn Number so that the flew always falls 
into the free molecular flow regime. 

3. A nozzle design procedure is developed to generate the nozzle 
geometry based upon given inlet total conditions that will 
produce droplets of the desired size. ’With constant expansion 
rate, higher total temperature and lower expansion rate would 
give flatter nozzle shape. Although different designs (nozzle 
geometry and inlet total conditions) can all reach the same goal, 
a state-of-art choice has to take into account of the boundary 
layer effect, multi-dimensional effect, precision of fabrication, 
strength of nozzle material etc. into consideration. Temperature 
dependence and variable expansion rate design can be implemented 
iteratively with the current proposed design procedure. 



31 


3.1 Recommendations 


To accurately solve the problem involving condensing flow, 
further investigations on the following subjects are recommanded : 

A. surface tension of small liquid droplet or surface free 
energy of small solid crystal with 10 to 100 molecules; 

B. an accurate nucleation theory and experiments which can 
accurately verify the calculated result of the theory; 

C. the state of the condensate belcw sublimation or triple 
line (solid or liquid); 

D. an accurate theory for rate of crystallization; 

E. values of thermodynamic properties such as enthalpy, 
density, specific heat etc. at all range of temperature 
and pressure. 

It is also suggested that experiments of test nozzle be conducted 


with nozzle design procedure described in the previous chapter to 
confirm the numerical computation presented in this work. 
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APPENDIX A 


DERIVATION OF DEICH'S NUCLEATION THEORY 


The nucleation rate consists of two components which we will 
derive separately as follows. 


1. First component 


This component is arised from the collision of many single 
molecules which happen to form the droplet of critical size. The 
number of molecules to make the critical size per unit mass and 
per unit time 1^ can be expressed as: 


I = - 



( A . 1 ) 


where (3 = number of molecules bombarding at a unit surface per 
unit time 

N = the equilibrium number distribution of the droplet 
z = f/N is a non-dimensional number distribution where f 
is a non -equilibrium number distribution 
n = number of molecules in each droplet 
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The expression for j3 is well known and can be found in any gas 
kinetic theory book to be 



The physical meaning of Equation (A.l) may be interpreted as 

follows: 4?tr 2 £ represents the number of molecules reaching the 

cr 

surface of a droplet of critical size per unit time. N is the 

number of droplets with critical size. The product of the two 

would be the total number of molecules impinging on all the 

droplets with critical size. In other words, this product 

represents the number of molecules condensing on the nuclei cf 

critical radius oer unit time, n is the number of molecules in 

cr 

a droplet with critical size. " ( ) c „ represents the change of 

the non-dimensional number distribution from the droplet with 


n - 1 molecules to the droplet with n molecules. Thus the 

cr cr 

product of these three terms gives the net molecule number flux 


growing from n - 1 molecules to n molecules, i.e. the 

cr cr 


nucleation rate. 
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According to Wegener [1.13] # equilibrium number distribution N 
can be expressed as 

2 

- / nkTA - 4rrr <r \ . * 

N = N exp ( — ) (A. 3) 

where the Gibb's enthalpy change associated with the formation of 

the droplet against the pressure difference (p - p ) to obtain 

the drop initially at p at the actual vapor pressure p was 

00 

neglected, and the ideal gas assumption was used. If we use R and 
N instead of m and k in Equation (A.. 2) and (A. 3) , also noting 
that the number of molecules n in a droplet with radius r can be 
expressed as : 

n = — 7 tr 3 p f N (A. 4) 


Equation (A. 2) and (A. 3) can then alternatively be expressed as: 


& = 


pN 


(A. 5) 


/ 


2ttRT 


4 3 — — 2 

-rrr p NRT A - 4trNr cr 

N = H exp ( ) 


(A. 6) 
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If we assume that 

1. The vapor has an equilibrium distribution before the 
expansion , 

2. The expansion is so rapid such that the number 
distribution of the droplet does not have enough time to 
adjust, thus remains at the initial equilibrium 
distribution, 

the non-equilibrium number distribution f would be the same as 

the initial equilibrium distribution N : 

0 

4 3 — — 2 

-?rr p NRT A - 47rNr <r 

f = N = N exp f — ) ( A . 7 ) 

O ' RT ' 

0 


With Equation (A. 6) and (A. 7), the non-dimensional number 
distribution z can now be expressed as: 


N 


L U f 

= — = — = exp ( 


4 3 — — 2 

RT A ^—7tr p N - 4 ttN r cr 
0 o 3 r t 0 

RT 


4 3 - - 2 

RTA-^ — irr p N - 4 tiN r <r 

— 


(A. 8) 


The derivative — can now be derived by chain rule, 
dn 
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dz _ dz dr 
on ” dr dn 


(A. 9) 


From Equation (A. 4), we have 


^ = 47tr 2 p N (A. 10) 

dr K f 


therefore, 


3r 

dn 


1 


2 - 
4nr p N 


(A. 11) 


Differentiate Equation (A. 8) with respect to r gives 


dz 

dr 


2 — — 2 — — 

RT A 47ir p N - 8TrNrcr RTA4rrr p N - 37rNrcr 

f 0 0 o f ^ 

l RT RT / 


exp ( 


. 4 3 - , - 2 

RT A — rrr p N - 4nNr <r 

0 0 3 0 

RT 


4 3 — - 2 

RTA_-7rr - 4nNr <r 

RT 


(A. 12) 


Now, substitute Equation (A. 5), (A. 6), (A. 9), (A. 11) and (A. 12) 


into Equation (A.l) we get 
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I = 
1 


pN 


/ 2 "* T ' p f 


RT A 4rcr~ p N - StiNt cr 

0 0 cr f cr 0 

RT 


7 _ _ 4 3 - „ — 2 

RTA4mr~ p N - 8nNr cr RT A — irr p N - 4nNr cr 

r r cr ^ / 0 0 3 cr f cr 0 

j exp [ 


RT 


RT 


4 3 - - 2 

RTA— r~irr p N - 4nNr <r 

3 cr f cr 

RT 


(A. 13} 


The number of nuclei with critical size generated per unit time 
per unit mass of condensing vapor can be computed by: 


I 

J = — — (A. 14) 

1 n 

cr 

where 

n = -4“ rcr 2 P N (A. 15) 

cr 3 cr t 


Substituting Equation (A. 13) and (A. 15) into Equation (A. 14) and 

with r expressed as 
cr 


2cr 

r cr ” p f RTA 


(A. 16) 


we can now obtain the final expression of the first component 


nucleation rate as 
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J 

i 


3 ANpv^RT 
g 

VSn o' 


TAir /cr - T A ] 


Z” 3TA<r /cr - 2(TA) 

g 0 0 0 

exp 

1 0 0 

T 

2 T A 

0 

v. * 


A 0 


(A. 17) 


where 


167tN<t 3 

3p 2 ( RT ) 3 
t g 


2 . Second component 

From chapter 2 we know that the droplets v/ith mean radius smaller 
than critical radius is unstable and will evaporated provided 
that the vapor condition (e.g. temperature and pressure) remains 
unchanged. However, during the expansion, the vapor condition 
will change and so dees the critical radius. We have shown that 
the critical radius is decreasing during the nucleation zone. 
Thus some of the droplets with subcritical size formed upstream 
would not have enough time to evaporate when they are transported 
to a downstream region, where they become supercritical. 
Consequently, these droplets would then become stable and serve 
as the nuclei for further condensation. This constitutes the 
second component of Deich's nucleation theory, which is not taken 
into account in the classical and Lothe-Pound nucleation 


theories . 
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It is obvious that the higher the expansion rate is, the faster 
will the rate cf decrease of the critical radius be. And the 
faster the rate of decrease of the critical radius, the larger 
the number of droplets would become supercritical fiom 
subcritical. Also changing of critical radius r^^ is the same as 
changing the number of molecules in a critical droplet n cr * 
Therefore, if the same two assumptions used in deriving the first 
component nucleation rate are made, as well, the number of 
molecules that v/ould make critical size droplet per unit time per 
unit mass of the condensing vapor I can be computed as: 


I 

2 


N 


0 


dn 

cr 

dt 


(A. 18) 


Using chain rule 1^ can be 


further expressed as: 


I 

2 


N 


0 


dn 

cr 


dr 

cr 


dr 

cr 

dAT * 


dAT 

dt 


(A. 19) 


Differentiate Equation (A. 15) we get 


dn 

cr 

dr 


4rcr 


p N 
K f 


cr 


(A. 20) 
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In addition to Equation (A. 16), can also be expressed as 


[1.06] : 


r 

cr 



(A. 21) 


With AT = T - T , r would look like 
s g cr 


r 

cr 


2cr 


P h 

K f fg 


M i ♦ r-) 


(A. 22) 


By differentiating Equation (A. 22) with respect to AT we would 
have 


dr 

cr 

dAT" 



(A. 23) 


Substitute Equation (A. 7), (A. 20) and (A. 23) into Equation (A. 19) 

and utilize the expression of r^ with Equation (A. 16) will give 


I = 


„ 3 -2 

4rcr p N 
c r f 


VC T 1 ) 


dAT 

dt 


exp ( i wW 


- 2 

4nNr o* 
cr 0 


RT 




(A. 24) 
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Similar to Equation (A. 14), 


can be expressed as 


I 

J = —L (A. 25) 

2 n 

cr 


By substituting Equation (A. 16) into Equation (A. 15) and (A. 24), 
the second component of Deich's nucleation theory can be 
computed as: 



The total nucleation rate J is the sum of the two components 


J = J + J 
1 2 


(A. 27) 



APPENDIX B 


RELATION BETWEEN SUPERCOOLING RATE dAT/dt AND EXPANSION RATE p 


As shown in Appendix A, the second component of the Dench's 
nucleation rate depends on the supercooling rate dAT/dt, which 
can be shown to be related to the expansion rate as follows. 


By definition AT = T - T , thus 


dT dT 

£AT _ £ £ (B . 1 ) 

dt dt dt 


Assume isentropic expansion, we have 


P = 



constant 


1 dp -y 1 9 

p dt % - 1 T dt 

^ g 


<B.2) 


(B . 3 ) 


thus 



( B . 4 ) 
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Also from Claus ius -Clapeyron equation we have 


where v is the 
sg 

p. With Equation 


, h 

ap :g 


aT v T 


(B.5) 


s sg s 

specific volume of saturated vapor at pressure 
(3.3) , this gives 


dT 

s 

dt 


v T 
sg s dp 


h . dt 


pv T 
sg s 



(B . 6 ) 


Hence, substituting Equation (B.4) and (B.6) back into Equation 
(3.1) we have 


pv T 

f T 1 1 1 )p 

dt l s r K } 

i-3 


(B .7 ) 


which shows that the supercooling rate dAT/dt is directly 
proportional to the expansion rate p. 



APPENDIX C 


DERIVATIVE OF TOTAL WETNESS 

To solve the equation of velocity, the expression of the 
derivative of total wetness dw/dx is required. Two ways we can 
derive this equation, mathematically or physically. Both methods 
give exactly the same expression. 

1. Mathematically 

dw/x can be obtained by directly differentiating w Equation (4.26) 


dw 

dx 


(x) 



W . ( x . ; 
1 i 


X) 


dx 

i 


(C.l) 


Assume 



i 


(C . 2 ) 


Substitute Equation (C.2) into (C.l), we have 
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ai a 


10 0 


w 

v 


( v ) 


d_ 

dx 


w. (x.; 

1 i 

x. =0 


x) dx 

i 


d_ 

dx 


_x 

x =0 



i 


x ) dx, 


d_ 

dx 


g ( x ; x ) 


g (0 ;x) 


6 g . , dx dg . , dx dg , . . dO dg 

= VMX;X — + -=2. x;x — - -=2- 0;x) -r- + ^( 0 ;x) 

dx dx dx dx dx dx ox 

Li j u i 


f^-(x;x) + ~{x;x) - |^(0;x) 

ox cx ox 


w { x ; x ) + ^2(x;x) - ^-(0;x) 
i ox ox 


Also from Equation (C.2) we know that 


^i. . , _ 3 _ \dg_ 

ox X i' X ox 5 x 


( x . ; x ) 


1 a , ,1 

■ sr la5 <3C i ; xl 

J i L 


Thus 


r 


x =0 
i 


5 w 

-5— -(x.; x) dx. 
5xi 1 


x =0 


5 

dx 


- *>] ^ 


dx 

dx 


(C. 3 ) 


(C. 4 ) 
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g(x.-xl - g(0)X> 


From Equation (C.5) and (C.3), we have 


dw , , 

d f 

3- (x) 

dx 

= d * Jx =0 


w ( x ; x ) dx 
i 1 i 


= w (x;x) + 
1 


3w , 

■=— (x.; x) dx, 

J„ 3x i i 


’x z0 


Differentiate Equation (4.24) with respect to 


Equation (4.2), we have 


dw 

-r-(x) 

dx 


— Xr | (X) P J (X) 

c(x) 3 cr f 


4tt 

c (x) 


r 2 (x ; x) a (x. ; x) J(x.) [T (x ; 

x i ii i ^ ^ 1 


:< =0 
i 


[h (T (x) ) 
g 9 


h (T (x ; x) ) ] 
f f , i i 


2. Physically 


(C.5) 


(C . 6 ) 


x and utilize 


x) - T (x) ] dx . 

g i 

c(x. ) 


(5.3) 


The increase of the total wetness can be attributed to two parts 
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(1) the birth of a new group of droplet, and (2) the growth of all 
the groups of the old droplets. The mass of the new born droplets 
per unit mass of the mixture is 

(dw) = p ~Trr 3 (x) J(x) dt (C.7) 

new f 3 cr 

The increase of mass of the group born at x^ currently located at 
x will be 

(dw) = p 4nr~(x ? x) dr J (x ) dt (C.3) 

old t i i i i i 


The increase of total wetness can then be calculated by 


dw 


(dw) + 
new 



i 


(C . 9 ) 


Mote that 


dt 


dx 
c (x) 


(C.10) 


Thus we have 


dw 


— (x) = p_ -nr 
dx e 3 cr 


(X) 


J (x) 
c ( X ) 


r 


dr J ( x ) 

i i 

p 4nr~(x ; x) dx. 

i i dx c(x.) i 

= 0 1 
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(C. 11) 


Also from Equation (4.2) we have 


dr 


p 4rtr (x ; x) 3 — (x. ; 
r z i i dx 1 


4jrr (x ; x) a (x ? x)[T. (x. ; x) - T (x) ] 

i i i i r,i 1 g 

[h (T (x) ) - h (T _ .(x.; x) ) ] c(x) 
g g f t , 1 1 


(C . 12 ) 


Hence, substitute Equation (C.12) into (C.ll) we obtain 


dw 

-r-(x 

dx 


1 _ 

C (x) 


4 

3 


Trr 3 (x) 
c r 


p c J(x) 

i. 


4tt 

+ C (x) 


r x 


X 


r 2 (x ; x) a (x ; x) J(x ) [T. (x. ; x) - T (x)]dx. 

i i 1 i 1 £ . 1 1 2 1 

(h (T (x) ) - h (T. . (x. ; x) ) ] c(x ) 

= 0 gg nr , 11 i 


(5.3) 



APPENDIX D 


DERIVATION OF EQUATION OF VELOCITY 

The equation of velocity can be derived from the differential 
form of continuity equation, momentum equation, energy equation, 
and the equation of state. These equations are listed as follows. 

continuity equation 

+ + + 1— = o ( d . l ) 

p c A 1 - w 


momentum equation 


cdc + 



dp = 0 


(D.2) 


energy equation 

/ 

cdc + (1 - v) c dT + wc dT = h dw 
p g f f fg 


(D.3) 


equation of state 

p = pRT 

g 


(D.4) 


First, differentiate equation of state Equation (D.4) gives 
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dT = 1 (i dp - — RT ) 
g R p P 9 


(D.5) 


1 dp 

Use Equation (D.2) for - dp , and Equation (D.l) for — in 

Equation (D.5) we get 


dT = - 

g R 


-cdc 


1 - w 


dc dA 
c A" 


r ac 

* ( - * T- • 


r -2S- 1 st] 

1 - w J gj 


(D.5) 


Substitute Equation (D.6) into Equation (D.3) gives 


cdc + (1 - w)c — 
o R 


-cac 

T 


+ (*£ + + ^L) RT ] 

- w v c A 1 - w' gj 


+ wc dT. = h _ dv 
t n z g 


(D.7) 


Group all the terms with dc together and divide the whole 

equation with (1 - w)c.T we have 

p g 


dc 

c 


re" (1 - c /R) 

o 


(1 - w)c T 
p g 


- l 


, wc dT - h dv 

^ + d 'C + — £ £ ( D . 8 ) 

A 1 - w (1 - w)c T 

p g 


Note that for ideal gas R = c - c , thus Equation (D.8) will be 

p v 


dc 

c 


(1 - w)?RT 


- 1 


wc dT - h dv 
dA dw t z f g 

A + 1 - w + (1 - w) c T 

p g 


(D.9) 
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From the definition of Mach number 




2 


c 

3rRT 


( D . 1 0 ) 


we can now obtain the equation of velocity 


dc 

dx 



dA 1 
dx A 


wc dT 

t i 

(l-w)c T dx~ 
p, g g 


c T -h # . 


p . g g 
(l-w)c 


T 

? 5 


dw 

dx 


( 5 . 2 ) 


Utilizing the differential form of the two component continuity 
equation, momentum equation, energy equation, and the equation of 
state, with exactly the same procedure, we can also yield the 
equation of velocity for two component model. The arithmetic for 
such derivation is straight forward and lengthy, thus will not be 


presented here. 



APPENDIX E 


DERIVATION OF MEAN MOLECULAR WEIGHT AND DENSITY 


1. Mean Molecular Weight 


The initial mass fraction of the condensing vapor w is defined 
as 


m m 

vo v: 

0 m m ^ + m 


(E.l) 


m = initial mass of condensing vapor 

m = mass of noncondensing inert gas 
n 

m = total mass of the mixture 


wetness w, the mass fraction of the condensate, can be evaluated 
as 


m m - m 

c vO v 

m m 


(E.2) 


where 

m = mass of condensate 

c 

m = remaining mass of condensing vapor after 

V 

ccndensat ion 
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Since the number of moles of the mixture vapor is the sum of the 
moles of the condensing vapor and inert gas, the mean molecular 
weight for the vapor phase defined as the mass per mole of the 
mixture vapor can then be related to the molecular weight of the 
condensing vapor M and inert gas M by 


m + m 
v n 


M 

g 


m m 
v n 

sr + m~ 

v n 


( E . 3 ) 


From Equation (E.l) we have 


m 

n 


m 


1 - w 

0 


( E . 4 ) 


From Equation (E.l) and (E.2) we have 


m 

— = w - w ( E . 5 ) 

m 0 


M can then be calculated from Equation (E.3) as 
g 


M 

g 


m + m 

v n 

m m 
v n 

m“ + m~ 

v n 


w - w + 1 - w 

0 0 


(w 


“> 5- * 11 

V 


V 


1_ 

M 
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1 - w 

(W - w) + (1 
0 M 

v 



(E.6) 


2 . Vapor Density 


Assume 

both 

condensing vapor 

ideal 

gases 

and of the same 

[A. 01] 




and noncondensing inert gas are 
temperature. From Dalton's rule 


P = P + P 

v : 


p — T 
H v M g 

V 


n M 



(E.7) 


Thus we obtain 


where M is 
g 


M 


P P 

v n 

M~ + m" 

v n 


expressed as Equation (E.6). 


(E.3) 



APPENDIX F 


INPUT OF THE FORTRAN PROGRAM 


1. Thermodynamic Properties 

PS AT (T) : saturated pressure P (T) as a function of 

sat 

temperature T 

TSAT(P) : saturated temoerature T (P) as a function of 

sat 

pressure P 

HG (T) : vapor enthalpy h (T) as a function of temperature T 

g 

TVAP(HG) : vapor temperature T (h ) as a function of vapor 

5 9 

enthalpy HG 

CP ( T) : constant pressure specific heat of vapor c^{T) as a 

function of temperature T, in some cases remains 

constant in the temperature range of interest 

GMISEN : ratio of specific heat of vapor c / c , in some 

p v 

cases it remains constant in the temperature range of 
interest 

HF (T) : liquid enthalpy h (T) as a function of temperature T 

SF (T) : liquid entropy s f ( T f ) as a function of temperature T 

DF(T,P) : liquid density p (T, P) as a function of temperature 
T and pressure P, in some cases remains constant 


110 



Ill 

in the temperature ar.d pressure range of interest 
CF(T) : specific heat cf liquid c^(T) as a function of 

temperature T, in some cases c ^ remains constant in 
the temperature range cf interest 
SIGMAO(T) : suface tension of flat surface cr^(T) expressed in 
the form of fitting polynomial NTERMS, SB(I), 

SC ( I ) , SD ( I ) 

MV : molecular weight M of the expanding vapor 


2. Nozzle Geometry 

AR ( x ) or PR ( x) : area ratio or isentropic pressure ratio as a 
function of st Teamwise distance of nozzle 
expressed in the form of coefficients of 
fitting polynomials NTERMA, AB(I), AC ( I ) , 

AD ( I ) or NTERMP, PB ( I ) , PC ( I ) , PD ( I ) 

N : No. of grids 

DIST : overall distance 

XO : starting location of the nozzle to the throat 


3 . Constants 
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R : universal gas constants 
NA : Avagadra Number 
PI : n 


4. Initial Conditions 

TO : inlet total temperature 
?0 : inlet total pressure 
HGV : inlet total enthalpy h. 
SOV : inlet total entropy s 

0 


5. Iteration Parameters 

NITER : maximum times of iteration 

RPSM, EPS : convergence criteria of the relative deviation 

of the iteration of Mach No., and other variables 
such as velocity VEL, vapor temperature TVAP , 
pressure P, droplet radius RI 
RELAX, RELAX 2 : relaxation factor of the intermediate 

iteration result before and after the end of 


nucleat ion 
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6. Experimental Data 

NEXP : No. of experimental data 
PREXP : experimental pressure ratio 

UNITS USED IN THE PROGRAM 


1 . 

Pressure 

: N / itT 



2 . 

Temperature : K 



3 . 

Density : 

Kg/m' 



4 . 

Velocity 

: m/s 



5. 

Length : 

m 



6. 

Enthalpy 

: J/Kg 



7 . 

Entropy : 

J/ (Kg K) 



8 . 

Specif ic 

Heat : J/ (Kg K) 



9 . 

Universal Gas Constant : 

J/ (Kg 

mol 

10 

. Molecular Weight : Kg/ (Kg mol) 


11 

. Surface 

Tension : N/m 



12 

. Thermal 

Conductivity : 

J/ (s m 

K) 



APPENDIX G 


DEFINITION OF THE CORRELATION COEFFICIENT 

The correlation coefficient can be defined as follows. Given n 

pairs of observations {x ,y ), the correlation coefficient r can 

i i 

be shown to be [G.01]: 



r serves as a measure of the degree of linear relationship among 
variables. The value of r is between 1 and -1. The larger the |r| 
is, the stronger the linear correlation is. 
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Ta b 1 e 6.1 Summa ry c f one c cir.pc ne r. - t e s t c a s e res v_i 1 1 
C: classical theory, D: Deich's theory 

error (Pa): standard deviation cf pressure ccmpariscn between 
theory and extrer iment 

STD ( 1 0 ~ ’ jim ) : standard deviation of droplet radius distribution 


Te s t Case 

Nucl 

O'/ O' 

00 

V K) 

f (Pa) 

?_ (Pa) 

error 

r ( Jim) 1 

:ti 

STD 

Barschdcr f f 

D 

.305 

376.9 

73323 

33246 

.0112 

.0241 

3.65 

No. 4 

C 

1.020 



33447 

.0124 

.0262 

4.29 

Barschdor f f 

D 

.310 

333.4 

73323 

27367 

.0172 

.0156 

3.03 

No. 7 

C 

.595 



27650 

.0179 

.0163 

3.51 

Barschdor f f 

D 

.340 

399.3 

73323 

21941 

.0023 

.0091 

2.4 5 

No . 1 0 

C 

1.000 



22036 

.0021 

.0093 

3 .39 

Binnie and 

D 

.365 

411.9 
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47309 
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.0251 

3.95 
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C 
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D 1 

.350 
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1.030 


! 

43711 

.0095 
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4.74 
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mm 
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1.44 
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i 



.0044 

1.56 

Gyarmathv and 

D 

.375 

420.7 

206340 

72633 

.0062 

.0634 

7.30 

Meyer No. 102 

C 

1.060 



72392 

.0053 

.0702 

3.54 

Gyarmathy and 

D 

.830 

392.3 

62052 

13346 

.0027 

.0154 

3 .30 

Meyer No. 108 

C 

.995 



18511 

. 0023 

.0164 

3.70 

Gyarmathy and 

D 

.805 

339.8 

62052 

13148 

.0073 

. 0057 
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C 
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.0070 

.0063 

1.96 
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D 
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333.7 
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C 

.905 
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mi 
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D 

.960 

442.1 

452230 

137080 

.0077 

.0636 

7 . 07 

Meyer No. 119 

C 

1.150 
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.0075 

.0740 

7 .36 
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D 

.330 


212010 

47295 

.0044 


2.62 

Meyer No. 91 

C 

.990 



47330 

.0046 


3.29 

Kadambi No. 3 

D 

.780 

352.3 

33503 

13664 

.0283 

.0066 

1.96 


C 

.995 



13742 

. 0290 

.0072 

2.42 

Kadambi No. 4 

D 

.750 

363.2 

25920 

14550 

.0298 

.0060 

2.20 


C 

.955 



14327 

.0291 

.0063 

2.32 

Kadambi No. 5 

D 

.785 


52332 

21295 

.0495 

.0036 

3.36 


C 

.935 



21565 

.0500 

.0091 

3.90 

Moore No. A 

D 

.735 

354.6 

24407 

9297 

.0228 

.0079 

3.41 


C 

.925 



9316 

.0231 

.0094 

2.57 
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TABLE 7.1 mean radius calculated at nozzle exit 13 cm downstream 
to the nozzle thcat for nozzles of various design 

STD ( 1 0 ~ 3 fim) : standard deviation of droplet radius distribution 


nozzle 

5 

O 

E-* 

p (1/sec) 

P (MPa) 

0 

r (/im) 
m 

STD 

NZl 1 

500 

2500 

1.7 

0.2795 

5.96 

NZ11 

500 

2500 

0.87 

0.0715 

12.3 

NZ12 

500 

5000 

1.7 

0.2117 

4.35 

NZ12 

500 

5000 

0.87 

0.0396 

6.23 

NZl 3 

50 0 

10000 

1.7 

0.1304 

4.33 

NZl 3 

500 

10000 

0 . 37 

0.0534 

5.50 

NZ21 

425 

2500 

0.33 

0.1339 

9.71 

NZ2 1 

425 

2500 

0.17 

0 . 0407 

Bfl 

NZ22 

425 

5000 

0.33 

0.0369 


NZ22 

425 

5000 

0 . 17 

0.0250 

3.86 

NZ23 

425 

10000 

0.33 

0.0451 

5.93 

NZ23 

425 

10000 

0.17 

0.0152 

2.37 

NZ3 1 

350 

2500 

0.027 

0.0111 

2.54 

NZ3 1 

350 

2500 

0 . 013 

0.0037 

1.22 

NZ3 2 

350 

5000 

0.027 

0.0073 

1.92 

NZ3 2 

350 

5000 

0.013 

0.0032 

.851 

NZ3 3 

350 

10000 

0.027 

0.0057 

1.00 

NZ33 

350 

10000 

0.013 

0.0026 

. 626 
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Figure 6.3 Droplet Size Distribution Of 
Barsohdorff No. 4 
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Figure 6.12 Droplet Size Distribution Of 
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Figure 6.18 Droplet Size Distribution Of 
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Figure 6.21 Droplet Size Distribution O t 
Gyermathy And Meyer No. 112 

Deich's classical 

6.20 Pressure Ratio Of Qyarmethy theory theory 

Meyer No. 102 For Deioh** Theory 
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Figure 6.24 Droplet Size Distribution Of 
Qyermethy And Meyer No. 108 
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Figure 6.2 7 Droplet Size Distribution Of 
Gyarmathy And Meyer No. 1 07 
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Figure 6.3 3 Droplet Size Distribution Ot 
Gyarmathy And Meyer No. 119 
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Figure 6.36 Droplet 8ize Distribution Of 
Qyermathy And Moyer No, 61 
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Figure 6.3 9 Droplet Size Dietribution Of 
Kadambi No. 3 
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Figure 6.42 Droplet 8*ze Dlatrlbutlo Of 
Kadambl No. 4 
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Figure 6.45 Droplet Size Distribution Of 
Kadambi No. 5 
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Figure 6.60 Nucl. Flux And Supercooling 
rtf Bsraohdorff No. 7 For Deloh’e Theory 
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Figure 6.52 Maoh No. And Entropy Ot 
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Fiflire 6.64 Velocity And Expansion Rato Figure 6.65 Knudaen No. And Prandtl No. 

Of Barachdorff No. 7 For Deloh'a Theory n-,.„hHnrff No. 7 For Delch'a Theory 
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Figure 6.86 Two Component* Of Nuol. Flux F j ouro 6 . 57 surface Tension Ratio Versus 

Of Baraohdorff No. 7 For Deloh’a Theory Total Temperature For Deich's Theory 
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Figure 6.63 Surface Tension Ratio Versus 
Total Entropy For Deich's Theory 
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Figure 6.64 8urf*ce Tenelon Retio Vereue 
Total Entropy For Cleaeloal Theory 
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Figure 6.7 3 Effect Of Total Pressure 

Figure 6.7 2 Effect Of Total Pressure On Droplet Size Distribution (Deioh's) 

Mean Radius (Deioh's Theory) __ 
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Figure 6.88 Two Component* Of Nucl. Flux 
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Figure 6.9 2 Velocity end Expenelon Rate Figure 6.9 3 Velocity And Expansion Rato 

Of Bareohdorff No. 7 (Two Comp., wO-t.O) of Barachdorff No. 7, Two Comp. (w0-0.6) 
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Figure 6.04 Velocity And Expansion Rets 
Of Barsohdorff No. 7 (Two Comp., wO-O.2) 
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